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Abstract 



Superintegrable systems are classical and quantum Hamiltonian systems 
which enjoy much symmetry and structure that permit their solubility 
via analytic and even, algebraic means. They include such well-known 
and important models as the Kepler potential, Calogero-Moser model, 
and harmonic oscillator, as well as its integrable perturbations, for exam- 
ple, the Smorodinsky-Winternitz (SW) potential. Normally, the problem 
of classification of superintegrable systems is approached by considering 
associated algebraic, or geometric structures. To this end, we invoke the 
invariant theory of Killing tensors (ITKT). Through the ITKT, and in 
particular, the recursive version of the Cartan method of moving frames 
to derive joint invariants, we are able to intrinsically characterise and in- 
terpret the arbitrary parameters appearing in the general form of the SW 
superintegrable potential. Specifically, we determine, using joint invari- 
ants, that the more general the geometric structure associated with the 
SW potential is, the fewer arbitrary parameters it admits. Additionally, 
we classify the multi-separability of a recently discovered superintegrable 
system which generalizes the SW potential and is dependent on an addi- 
tional arbitrary parameter k, known as the Tremblay-Turbiner-Winternitz 
(TTW) system. We provide a proof that only for the case k — ±1 does 
the general TTW system admit orthogonal separation of variables with 
respect to both Cartesian and polar coordinates. 
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Chapter 1 



Introduction 



The origins of this thesis can be traced back to 1965 when Winternitz and Fris 
used the theory of Lie groups to classify the separable webs of the Euclidean plane. 
Specifically, they computed the invariants of the second-order symmetries of the 
Laplace equation 



d 2 V | d 2 V 
dx 2 dy 1 

defined on the Euclidean plane under the action of the Euclidean group. As we 
shall develop in Section 2.3, the classification of such separable webs is intimately 
related to the Hamilton- Jacobi theory, which provides a formulation of Hamiltonian 
mechanics that allows one to conclude the solubility of a classical Hamiltonian system 
(see Section 2.1). Indeed, by determining the orthogonal coordinate systems that 
afford separation of variables in the corresponding Hamilton- Jacobi equation for a 
Hamiltonian system (see Section 2.2), one can ultimately obtain exact solutions, or 
trajectories, of the Hamiltonian system. 

The work of Stackel in 1893 and its extension by Eisenhart in 1934 were both 
essential to the development of this geometric approach to studying Hamiltonian 
systems. Both of these canonical papers established the essential role of a geometric 
object that naturally arose when studying orthogonal separation of variables of the 
Hamilton- Jacobi equation: a Killing tensor of valence two. Studying the intrinsic 
properties of Killing tensors corresponding to certain orthogonally separable cases 
thus provided a new approach to the problem of solving a Hamiltonian system. In 
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2002, McLenaghan, Smirnov, and The [45] utilised classical invariant theory to de- 
velop the properties of such Killing tensors, which marked the initial stages in the 
development of the invariant theory of Killing tensors. Since then the theory has 
been utilised by a significant number of authors who have applied it in a more gen- 
eral setting to study Killing tensors defined in spaces of constant curvature (e.g., see 
[45, 17, 1, 3, 13, 12, 15, 29, 30, 31, 32, 55, 56, 70, 73] and the relevant references 
therein) . 

In this thesis we use the invariant theory of Killing tensors in the context of solving 
the classification problem for superintegrable Hamiltonian systems defined on the Eu- 
clidean plane 1 . One system of particular interest is the Tremblay-Turbiner-Winternitz 
(TTW) system introduced in 2009 [65], which is a generalisation of another super- 
integrable system of interest, the Smorodinsky-Winternitz (SW) system. The unique 
feature of the TTW system is that it provides an infinite family of solvable and in- 
tegrable (quantum) systems on the Euclidean plane with respect to a parameter k 
[65, 66, 67]. Furthermore, this family of systems was proven to be both classically 
and quantum superintegrable for any rational value of k [33, 35, 36, 37]. 

A Hamiltonian system with n degrees of freedom is completely integrable if in 
addition to its Hamiltonian, H, it admits n — 1 first integrals of motion, Fi,i = 
1, ... ,n — 1, that are well-defined functions on the phase space. These first integrals 
must be functionally independent, 

dH A dFi A ... A dF„_i ^ 0, 

and in involution with respect to the Poisson bracket (see Section 1.3) 

{H, F,} Po = {F u F^p, = 0, i,j = 1, . . . , n - 1. 

If such a completely integrable system admits additional n — 1 functionally indepen- 
dent first integrals of motion, Gi, i = 1, . . . n — 1, namely 

dH A dFx A ... A dF n _x A dC?i A ... A dG n _i ^ 0, 
1 The term superintegrable was introduced in 1983 by S. Wojciechowski [71]. 
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which are in involution 

{H,Gi} Po = {Gi,Gj} Po = 0, i,j = 1, .. . ,n - 1. 

then the system is maximally superintegrable (see Section 3.1). Since such systems 
have more constants of motion than degrees of freedom it is possible to integrate them 
through analytic and algebraic methods. For this reason, superintegrable systems are 
often sought after as starting points for studying much more complicated systems. 
Some familiar examples of superintegrable systems (in the realm of physics) include 
the harmonic oscillator, the Kepler system, the Calogero-Moser system, the hydrogen 
atom, and others (see e.g. [7, 21, 46, 48, 58, 62]). 

The construction of an infinite family of superintegrable systems, such as the TTW 
system, is fairly straightforward. We demonstrate the technique with respect to the 
TTW system, which begins with the second-order system defined by a Hamiltonian 
(see Section 2.1), of the form 

H = pl+p 2 y -^ + y*) + ^ + l (1.1) 

where p x ,p y denote the momenta coordinates, and x, y denote the position coordi- 
nates. This is the SW system, which is recognised as a seond-order perturbation of 
the two-dimensional harmonic oscillator and further, is an example of a superinte- 
grable Hamiltonian system since it is multi- separable with respect to both Cartesian 
and polar coordinates (see Section 2.3). The property of being superintegrable is 
what we intend to retain in the construction of a new Hamiltonian system. 

To this end, we can transform the Hamiltonian (1.1) to polar coordinates and 
replace 8 with k9, where k is some new arbitrary parameter. Such a transformation 
of course preserves the property of the system being separable with respect to polar 
coordinates, though we now arrive at a new Hamiltonian of the form 

*=rf+^ + ^--V + I(-^ + -| s ). (1.2) 

This is the TTW system. One characterising feature of the TTW system is that much 
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like the system (1.1), it is still separable with respect to polar coordinates and it still 
behaves as a perturbation of the harmonic oscillator. Though this development raises 
the question of whether the system remains superintegrable for all values of k. In their 
initial publication, Tremblay et al. [65] conjectured and provided strong evidence for 
the superintegrability of this system for all rational values of k, a claim which fueled 
the activity of researchers who were set out to prove its validity. Indeed, if the system 
remained superintegrable then it provided an example of a second order system that 
could generate an infinite number of higher order superintegrable systems, which 
would work towards the development of a classification theory for superintegrable 
systems. The claim was proven to be true in the classical case, where Kalnins et 
al. [37] explicitly solved for the higher order first integrals, appearing as polynomials 
in the momenta, by using the structure algebra of first integrals. In this context, 
one must determine the order defining when the algebra of first integrals closes with 
respect to the Poisson bracket, and find the corresponding structure equations for the 
symmetry algebra (see [34] for more details). It was also proven in the quantum case 
through a recurrence method introduced by Kalnins et al. (see e.g. [34, 35, 36]). In 
this thesis we investigate whether there are values of k for which the system admits 
two first integrals of motion quadratic in the momenta. Prior to this thesis, such a 
result has not been shown. 

In the Euclidean plane, maximal superintegrability requires that a Hamiltonian 
system defined by H (see Section 2.1) admits two additional first integrals of motion, 
Fi,i = 1,2 functionally independent and in involution, meaning that 

dH A dF 1 A dF 2 ^ 

and {H,Fi} = Xn{Fi) = 0, for i = 1,2. Note that a Hamiltonian is always a 
first integral of motion, a property which can be physically interpreted as a system 
demonstrating conservation of energy. For the system given by the Hamiltonian (1.2), 
we also readily obtain a second first integral by exploiting the fact that the system 
is separable with respect to polar coordinates. Indeed, since (1.2) takes the form of 
a natural Hamiltonian, we can use a result by Liouville [42] which yields that such a 
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system will admit a first integral of motion quadratic in the momenta according to 

F(q\ Pi ) = ^(fipiPj + U(q l ), i,j = 1,2, (1.3) 

with position coordinates q\ momenta coordinates pi, where is a Killing two- 
tensor whose normal eigenvectors generate the polar coordinate web (see Section 
2.3), and U(q l ) is a potential satisfying dU = KdV, where K = Kg -1 , i.e. K is a 
(1, 1)— tensor, which in component form is given by K % j := K %l gij. And so, in fact, 
there remains only one first integral of motion to find in order for us to conclude 
superintegrability of the TTW system. 

1.1 Summary of Results 

In this thesis we provide a new perspective on the study of joint invariants of Killing 
tensors, which are objects originally defined by Smirnov and Yue in 2004 [56], and 
further extended by Adlam et al. [3]. In particular, we establish a relationship be- 
tween the geometric and analytic properties of superintegrable Hamiltonian systems, 
defined on the Euclidean plane, E 2 , through the vanishing of joint invariants in the 
product space /C 2 (E 2 ) x /C 2 (E 2 ), where /C 2 (E 2 ) denotes the vector space of Killing 
two tensors defined on the Euclidean plane (see Section 2.4). We focus our study on 
the characterisation of the Smorodinsky-Winternitz potential, where we can readily 
establish a link between the arbitrary parameters appearing in the potential and the 
parameters in the associated vector space of Killing tensors. In particular, we find 
that the more structure we place on the product space /C 2 (E 2 ) x /C 2 (E 2 ), the more 
general form a corresponding superintegrable potential is allowed to take, and vice 
versa. 

As a generalisation of the Smorodinsky-Winternitz potential, we then focus our 
attention on the Tremblay-Turbiner-Winternitz potential, where we provide a defini- 
tive answer to the following question: 

For which values of k is the TTW system multi-separable? 

Through the geometric description of a superintegrable potential provided by certain 
types of Killing two-tensors (see Chapter 3.1), we find that only when k — ±1, which 
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in fact reduces the system to the Smorodinsky-Winternitz potential, does the TTW 
system retain the property of being multi-separable, in particular, with respect to 
both Cartesian and polar coordinates in the canonical position. 

1.2 Overview and Historical Development 

We begin with a Hamiltonian system defined (on a pseudo-Riemannian manifold A4) 
by a natural Hamiltonian of the form 

H{q\Pi) = y ij PiPj + V(q% i,j = 1,2, (1.4) 

where (q,p) denote generalised canonical coordinates. A first integral of a Hamil- 
tonian system is defined as a smooth function F = F(q,p) which is constant along 
the Hamiltonian flow. The flow can be interpreted as the geometric manifestation 
of Hamilton's equations in classical mechanics (see Section 2.1), and so one is often 
interested in the explicit form of a Hamiltonian's first integrals of motion. 

From a physical perspective, the first integrals characterise quantities that are 
conserved throughout the motion of a Hamiltonian system, a notion which evokes 
some familiar examples of constants of motion such as energy, and both angular and 
linear momentum. In a way, the more first integrals a system has, the more symmetry 
and structure the system must have. If we can determine n functionally independent 
first integrals in involution for a system with n degrees of freedom, then the system 
is said to be completely integrable, which means we can explicitly determine the 
trajectories or orbits of the system (see e.g. [62]). 

In general, an additional geometric structure is required to determine the first inte- 
grals of a Hamiltonian system. Some methods which have been successfully employed 
include the Lax representation method (see e.g. [51, 26]) and the bi-Hamiltonian ap- 
proach (see e.g. [53, 63] and references therein). In this thesis, we will make use of 
orthogonal separation of variables in the context of the Hamilton- Jacobi theory to 
identify quadratic first integrals of motion, and vice versa (see Section 2.3 for details). 
This approach to determining first integrals has a well-established history which we 
will discuss in Chapter 2. Furthermore, in Section 2.3 we will present the natural 
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link between orthogonal separation of variables and Killing two-tensors denned on 
the Euclidean plane. 




Figure 1.1: The orthogonal coordinate webs in the Euclidean plane, in the canonical 
position, arranged according to their number of singular points. Clockwise from top 
left: elliptic-hyperbolic, parabolic, polar, and Cartesian. 

Formally, the search for superintegrable systems in classical and quantum mechan- 
ics (in the context of this thesis) also began in 1965, when Fris et al. [22] published 
a list of superintegrable potentials defined on the Euclidean plane. Two years later, 
the study of superintegrable systems defined on three-dimensional Euclidean space, 
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E 3 , was initiated by Winternitz et al. [69] and Makarov et al. [43]. In 1990, these 
pioneering works received attention from Evans [21], who was able to extend their 
results by performing a systematic search for all maximally superintegrable potentials 
on E 3 that admitted (at most) quadratic first integrals. 

Evans claimed his list of second-order superintegrable systems to be exhaustive 
up to the equivalence class of linear transformations [21]. Though, there was one 
superintegrable system which did not find itself on Evans' list: the Calogero-Moser 
model. If we let (g 1 ,^ 2 ,^ 3 ) = (x,y,z) denote canonical Cartesian coordinates, then 
the Calogero-Moser system is defined in E 3 by a natural Hamiltonian, i.e. 

H= l(p 2 x +p 2 y +p 2 z ) + V(x,y,z) 

with potential given by 

V{x, y, z) = - — l —— + - — l —— + - — — — , (1.5) 
[x — y) 2 (y — z) 2 [z — x) 2 

which was well-known to be superintegrable at the time [71]. 

The reason such a model was overlooked by Evans' method was due to the fact 
that he assumed the associated Killing two-tensors would be in the canonical form 
[2]. This assumption in [21] and [22] excluded the possibility that a Hamiltonian 
system could be separable in an orthogonal coordinate system which was in a non- 
canonical position. We direct the reader to Figure 1.2 for examples of such coordinate 
systems on the Euclidean plane. Indeed, the Calogero-Moser system is an example 
of a system that admits five first integrals of motion of the form (1.3), where the 
associated Killing tensor that comes from solving the compatibility condition gives 
rise to five non-canonical characteristic Killing tensors (see [31, 2] for more details on 
the Calogero-Moser system; see Section 2.3 and 2.4 for the theory). 

Naturally, this motivates a generalisation of the aforementioned approach which 
makes use of non-canonical Killing tensors. In 2005, Adlam [1] began this endeavor 
and provided a more complete list of such superintegrable potentials defined on the 
Euclidean plane. Nevertheless, the case when a potential admits orthogonal separa- 
bility with respect to elliptic-hyperbolic coordinates in the canonical position, and 
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polar coordinates in a general position has not been completely classified 1 . 

As such, we seek to further develop the growing collection of superintegrable 
potentials defined on E 2 by providing a new interpretation of their functional form 
characterised via the invariant theory of Killing tensors. In particular, we derive 
a relationship between the joint invariants defined on the product space of Killing 
two tensors (consisting of Killing tensors compatible with a potential V), and the 
arbitrary parameters of a superintegrable potential which is multi-separable in two 
orthogonal system of coordinates. 

1.3 Notation and Conventions 

In this thesis, we let Ai denote an n-dimensional pseudo-Riemannian manifold, mean- 
ing that Ai is a smooth manifold endowed with a non-degenerate covariant symmetric 
metric tensor g. The tangent bundle of Ai, T(Ai), is defined as the disjoint union of 
the tangent spaces at all points of Ai: 

T{M) = [J T P (M), 

where the tangent space T p (Ai) is the set of all linear maps X : Ai — > R satisfying 
for all f,g E At 

X(fg) = f(p)Xg + g( P )Xf, 

i.e. the set of all derivations of Ai at p. An element of T p (Ai) is called a tangent 
vector. Whence, we also have the cotangent bundle of Ai, defined as the disjoint 
union of the cotangent spaces at all points of Ai: 

T*(M) = □ t;(m), 

where the cotangent space is the dual space of T p (Ai). Recall that a covector is 
a real- valued linear functional on a space, i.e. a linear map u : T p (Ai) — > M. A 

1 The choice of which coordinate system is in the canonical position is irrelevant, only that one 
system is in a general position. 
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Figure 1.2: Examples of orthogonal coordinate webs in the Euclidean plane, in some 
general position. Clockwise from top left: elliptic-hyperbolic, parabolic, polar, and 
Cartesian. 

significant fact which helps to identify T*(A4) ioi p E M is then that given a basis 
for the tangent space, (-Ej), the covectors (e J ) defined by 

J ' \ ifi^j 
form a basis for T*(A4), p G M. 
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Our results will take place in the Euclidean plane, so in general we will take 
A4 = E 2 where the metric tensor's components then become trivial, e.g. in Cartesian 
coordinates the components of the metric tensor are g^ = 5ij where 5ij denotes the 
familiar Kronecker delta. The remaining content of this section is devoted to the 
presentation of some specific definitions and theorems relevant to the later sections. 



1.3.1 Tensors and Index Notation 

Following the presentation in [40], a multilinear 1 function 

T : T*{M) x ... x T*{M) x T P {M) x . . . x TJM) ->• R 
^ v ' v v ' 

s copies r copies 

defined at a point p 6 M which maps r covectors and s vectors to a real number 
is called a tensor of valence (r, s), denoted TJ. In the case that r or s is identically 
zero, we will refer to such a tensor as a covariant s-tensor or contravariant r-tensor, 
respectively. As a few familiar examples, every covector uj : T p (Ai) — > K is identified 
as a covariant 1-tensor, the dot product on W 1 is a covariant 2-tensor, as a bilinear 
form, and the determinant, considered as a function on n vectors, is a covariant 
n-tensor on M. n . 

In Section 3.3 we will require all tensor components to be written in polar coordi- 
nates, and so we will make use of the usual tensor transformation law to acquire the 
appropriately defined objects. For example, the metric tensor defined with respect to 
a new set of coordinates (q 1 , . . . , q n ) is determined with respect to the old coordinates 
(q\...,q n ) via 

9ij = A fc jA jgw, 

where denotes the metric tensor in the new coordinates, and the A k i denotes the 
usual Jacobian matrix defined according to 



1 A map is said to be multilinear if it is linear with respect to each of its arguments, e.g. a 
multilinear function in two variables is a bilinear function. 
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relating "new" and "old" coordinates. 

As already noted above, we will frequently make use of the Einstein summation 
convention for tensor indices, by which any repeated upper and lower index implies 
summation over that index from 1 to n. As an example which will prove useful 
in Section 2.4, we have the components of a general Killing tensor defined on the 
Euclidean plane in Cartesian coordinates (g 1 ,? 2 ) = (x,y) as 

K ij = A H + 2e (^Y + Cei m eVV, 

where 

A«=( Pl M, B*=( M, C = K 

\ 03 02 J \ -05 J 

and e % j = g ie eej, where e^- denotes the familiar two-dimensional Levi-Civita symbol 1 . 
Indeed, making use of the Einstein summation convention we can explicitly write out 
the four components of this general Killing tensor as follows 

K 11 = A 11 + 2e\B\ 2 + CeV 2 gV 

= 0i + 2/% + /V, 
K u = A 12 + e\B 2 q 2 + e\BW + Ce 1 2 e 2 1 q 1 q 2 

= 03 - 05V ~ 04% ~ 06%y, 

K 22 = A 22 + 2e 2 1 B 2 q 1 + Ce\e 2 1 q 1 q 1 
= 02 + 2/3 5 x + /3 G x 2 . 



We have also made use of round brackets enclosing tensor indices to imply symmetri- 
sation according to 

= 2~(Tij + Tji), 

In an analogous manner, anti-symmetrisation is denoted by the use of square brackets 



-^k?1 2^ij Tji). 

In particular, this is a pseudotensor which is antisymmetric under the exchange of any two 
slots: eij = — e Jl . Its values are given by eia = \/\g\, where \g\ = det gij . [68] 
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Lastly we will discuss a few basic tensor operations. Tensors of the same valence can 
be added or subtracted by summing corresponding components in the anticipated 
way. Multiplication of tensors introduces the tensor product operation as follows. 
If T G Tg^{M) and S G TJ^(A4), then their tensor product T S is a new tensor in 
T S1 ^(M) defined by 

(T®S) (wi,... 

= T(wi, . . . , w ri ,Xi, . . . , X sl )S[T]i, . . . , r] r2 , Yi, . . . , F S2 ). 

Of particular importance in Section 2.4 will be the symmetric tensor product. If 
T\ G (Al) and T 2 G 7^ 2 (Al), then their symmetric tensor product is a symmetric 
tensor T\ T 2 given by 



where S n denotes the symmetric group on a set of n elements. It can be shown 
that the symmetric tensor product is associative, commutative, and bilinear. An 
analogous construction is defined by the wedge product A of tensors, which yields an 
antisymmetric tensor (see e.g. [40]). 

We remark that the lowering and raising of tensor indices will be done with 
the covariant metric and its inverse, the contravariant metric, respectively. For 
example, given a vector V with components V k relative to polar coordinates in the 
Euclidean plane, its corresponding 1-form V has components Vi given by 



(T 1 QT 2 )(X 1 ,...,X S1+S2 ) 




Vi = g ik V\ 



where g^ is the metric tensor in polar coordinates, i.e g 



( 



1 

r 2 



. In this case 



we see that V\ = V 1 and V 2 = r 2 V 2 . 
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1.3.2 The Lie Bracket and its Generalisation 

The covariant derivative defines how a vector field changes, i.e. it is a generalisation 
of the directional derivative from vector calculus. It can be defined as a connection 
V on the tangent bundle of a manifold, where a connection provides a prescription 
for moving vectors from one point to another on the manifold. The definition of a 
connection on a manifold does not require that a metric be defined on the manifold, 
but in the case that a non-degenerate metric is defined on a manifold, then there 
exists a unique torsion-free connection called the Levi-Civita connection. This choice 
of connection is defined by Vg = 0, i.e. the choice of connection with respect to which 
the covariant derivative of the metric vanishes. The components of the Levi-Civita 
connection are given by the Christoffel symbols, 



so that with this choice the covariant derivative of a covector field oji is given by 



where a comma followed by an index, e.g. denotes partial differentiation with 
respect to the coordinate qK The covariant derivative of a contravariant vector field 
V 1 is given by 



This can be generalised for tensor fields of higher valence (see e.g. [40]). If we work 
in Cartesian coordinates for E 2 then all components of the Levi-Civita connection 
vanish, and so all covariant derivatives reduce to simple partial derivatives. Accord- 
ingly, the covariant derivative gives us a way to "compare" vectors in neighboring 
tangent spaces, i.e. it gives a generalisation of parallel transport (see e.g. [40]). The 
Lie derivative is then another type of derivative between vector fields that one can 
define on a manifold. In particular, it is a canonical way of carrying out such a 
comparison that does not depend on the existence of a connection. Explicitly, the 





v' :; = Vj-v* = v\j + r jk v k . 
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1.3.2 The Lie Bracket and its Generalisation 



Lie derivative of a vector field Y along the flow of X is given by 

C X Y = X»Y V ^ - Y»X V ^ = [X, Y] u , 

where [ , ] denotes the Lie bracket which defines a new vector field as the commutator 
of two vector fields. It can be shown that the Lie bracket is bilinear, anti-symmetric, 
and satisfies the Jacobi identity. 

We now introduce the Schouten bracket and discuss a few of its properties. The 
Schouten bracket for contravariant tensor fields is a generalisation of the Lie bracket 
for vector fields. Formally, the Schouten bracket is defined as a real bilinear operator 
[ , ] : T p (Ai) x T q (M) T p+q ~ 1 whose operation on a pair of contravariant tensor 
fields P G T P (M), Q G T q (M) is defined by [54] 



[P Qy i --- i p+i- 1 = ^ y pih—ik-iWk—ip-iQ^Qhi—ip+q-i) 
k=l 

P 

_|_ ( — 1)^ pi i i--- i k-l\L 1 \ik---ip-iQ^Qip---ip+q-i] 

k=l 
1 

_ ^ ^ Q(ii...it_ 1 \[i\ii...i q -iQ^pi q ...i p+ q- 1 ) 

1=1 
1 

_ ^_iyq+p+q+£Q[h---k-i\v\k---iq-iQ^piq---i P +<i-i] 



where any index located between vertical bars means that the index is excluded from 
any anti-symmetrisation or symmetrisation. If we take P and Q to be tensors of 
type (p, 0) and (q, 0), respectively, then the Schouten bracket [54] of P and Q yields 
a tensor, [P, Q], of type (p + q — 1, 0). As an example, if we let P and Q be tensors 
with p — 1, so that P is now a vector field, and allow q to remain arbitrary, then by 
definition of the Schouten bracket we have that 

[P 1 Q} k ^ k " = (CpQ) kl - kq . 

And so indeed, when q — 1, we see that we obtain the Lie bracket of the vector fields 
P and Q. 
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1.3.3 Poisson Manifolds 



In general, a Killing tensor field of valence p defined on (Ai,g) is a symmetric 
(p, 0) (contravariant) tensor field K satisfying the generalised Killing tensor equation 
given by 

[g,K} = 0, (1.8) 

where [ , ] denotes the Schouten bracket. Since the Schouten bracket is M— bilinear, 
the set of solutions to the system of overdetermined PDEs given by (1.8) forms a 
vector space over R. We denote the vector space of valence p Killing tensor fields 
defined on M by JC P (M). 



1.3.3 Poisson Manifolds 

Consider the commutator of respective Hamiltonian flows in the classical Poisson 
bracket (see Section 2.1). First we define a Poisson bivector as a (2,0)-tensor, 

defined on M. which satisfies 

[Po,P o ] = 0, 

where [ , ] denotes the Schouten bracket. Any smooth manifold which admits a 
Poisson bivector, or the general Poisson bracket defined for smooth functions f,g : 
M ->■ K via 

U n\ - P« 9f 99 

is called a Poisson manifold. If we identify Pq as the canonical Poisson bivector, 
then the previous equation becomes 

P ° ^r* dpi dq i dpi dq i ' 

Following [40] , we can construct a smooth vector field on Ai for some smooth function 
H : T*(M) -)> R via 

X H = [P ,H], 
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1.3.3 Poisson Manifolds 



where the flow of this Hamiltonian vector field is called its Hamiltonian flow. With 
respect to (local) coordinates x a = (q\pi), the coordinates of the flow are then given 

as 

arr 

H dx a 

A Poisson manifold endowed with a smooth function H defines a Hamiltonian system, 
where the function H is aptly denoted as the Hamiltonian of the system. The integral 
curves of the Hamiltonian flow are then called the orbits or trajectories of the system. 

The Lie bracket and Poisson bracket can be shown to admit the following rela- 
tionship: 

Proposition 1.9 (cf. [40]). Let F,G : Ai — >■ R be smooth functions, where AA is a 
Poisson manifold. Denote their respective Hamiltonian flows by Xp and Xq- We 
then have that 

[X F ,X G ] = — X{ FjG }, 

where [ , ] is the Lie bracket from the previous section, and { , } is the Poisson 
bracket. 

A proof is provided in [40]. 
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Chapter 2 

Theoretical Overview 



In this chapter we provide the reader with the mathematical framework on which 
this study is built. This begins with a brief review of the Hamiltonian formulation 
of classical mechanics, which is the context of the problems we are interested in 
studying. We define and discuss what it means to solve a Hamiltonian system, which 
motivates our review of the Hamilton- Jacobi (HJ) theory of separability as a powerful 
method for solving particular Hamiltonian systems. The HJ theory will allow us to 
explicitly define what is meant in saying that a system is separable with respect to a 
system of coordinates on a smooth manifold, where we focus on results relevant to the 
Euclidean plane. In this context, we present theorems by Jacobi, Eisenhart, Liouville, 
and Benenti which establish the link between orthogonal separation of variables and 
Killing two-tensors defined on E 2 . This will motivate the final section, where we 
present the reader with a review of the invariant theory of Killing tensors. One goal 
in this final section is to introduce the method of moving frames by using a recursive 
version of it to derive the fundamental invariants of the vector space of Killing tensors 
of valence two under the action of the isometry group on E 2 , SE(2). 

2.1 Hamiltonian Mechanics 

In the Lagrangian formulation of classical mechanics the state of a mechanical system 
defined on an n-dimensional manifold Ai is completely determined once we specify its 
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2.1 Hamiltonian Mechanics 



generalised position and velocity coordinates. In the most general formulation, this 
is accomplished by defining a function, L(q t ,q l ), of the position coordinates q l and 
their velocities q % , where the dot denotes a derivative with respect to time, denoted 
by t. This function L is called the Lagrangian of the system and essentially defines 
a function from the tangent bundle on a manifold A4, denoted by T(M), to M. 

The evolution of a system can then be determined by employing Hamilton's prin- 
ciple or the principle of least action. We take the functional (see e.g. [39]) 



to define the action of the system, and from it compute the functional derivative 
of 5* with respect to each of the q % coordinates. By Hamilton's principle we have 
that the coordinates evolve in such a way (as functions of t) for which S remains 
stationary with respect to variations in q % that leave the initial and final time values 
unaffected. In other words, the action should satisfy bS = 0, where 6 represents 
the variation with respect to q % . Effecting this computation leads to the well-known 
Euler- Lagrange equations, or the equations of motion, for a given system, namely 



In the Hamiltonian formulation of classical mechanics, we consider a slightly dif- 
ferent set of generalised coordinates. In particular, we derive the canonical momenta 
coordinates pi from the Lagrangian L via 



In this way, the canonical momenta coordinates can be taken as defining a covector 
field in the cotangent bundle T*(A4). The relationship between the Lagrangian and 
Hamiltonian formulations is then given by the Legendre transformation between their 
respective coordinates, namely a smooth function H : T*(M) — > W defined by 




(2.1) 




(2.2) 




(2.4) 
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2.1 Hamiltonian Mechanics 



where, for the purpose of this thesis, we have assumed that H has no explicit time 
dependence. This function H is called the Hamiltonian of the system. Through 
equation (2.3) we see that the time derivatives appearing in (2.4) can be taken as 
functions of the generalised position and momenta coordinates. Indeed, substituting 

(2.3) into the Euler-Lagrange equations (2.2) and taking the partial derivative of 

(2.4) with respect to q l yields the evolution of our momenta coordinates, 

dH , . 

a = ~ W (2-5) 

Similarly, we can take a partial derivative of (2.4) with respect to p,i to obtain the 
evolution equation of our position coordinates, namely 

dpi 

Together, the equations (2.5) and (2.6) are called Hamilton's canonical equations. 
An important observation to make is that Hamilton's equations form a set of 2n 
first-order differential equations for the 2n unknown functions q l (t) and Pi(t) of the 
system. This system of first-order differential equations thus replaces the n second- 
order differential equations (2.2) admitted by the Lagrangian treatment. 

Hamilton's equations define a vector field on the cotangent bundle T*(Ai). The 
2n-dimensional phase space of a Hamiltonian, admitted by generalised position co- 
ordinates q % and generalised momenta coordinates pi, is an example of a symplectic 
manifold, which is defined as a smooth manifold with a smooth, non-degenerate, 
closed 2-form. In particular, we can construct such a 2-form via 



co = 

i=l 



Y^d Pi Adq\ (2.7) 



where d denotes the exterior derivative and A denotes the exterior, or wedge, product. 
This is the canonical symplectic 2-form. A fundamental result in the theory of sym- 
plectic structures, due to Darboux (see e.g. [40]), states that for any 2n-dimensional 
symplectic manifold with symplectic form u> there exist local canonical coordinates 
such that oj assumes the canonical form (2.7). Such coordinates are called Darboux 
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2.1 Hamiltonian Mechanics 



coordinates, or canonical coordinates. Further, the canonical symplectic 2-form u 
induces a pairing between the tangent and cotangent bundle, namely a bivector field 
on the manifold called the canonical Poisson bivector, which is obtained by taking 
the inverse of the canonical symplectic form 

Using the Poisson bivector, we can define the Poisson bracket for two smooth func- 
tions / and g defined on the cotangent bundle, by 



{f,9U = Pl 



.tj df_dg_ 
dq i dp 



3 

dj^dg_ _ dg df 
dpi dq % dpi dq i 



Given a symplectic manifold Ai with a canonical Poisson bivector Pq (2.8), if we have 
a smooth function H : T*(M) — > R, then 



X H = [P ,H] 

dh 

dpi dq % dq l dpi 



^dH d dH d . , , , , 
> ~ — 7— — — — r — — , (in local coordinates) (2.9) 

< J rln- Art 1 Art 1 rtn- 



defines the Hamiltonian vector field, where [ , ] denotes the Schouten bracket defined 
in Section 1.3. A symplectic manifold M. with a Hamiltonian vector field defined 
with respect to a Poisson bivector, and smooth function H : T*(A4) — > R, namely 
{A4, Xh, H}, comprise a classical Hamiltonian system. This development also yields 
the condition for a smooth function F : T*(A4) — > M to be an integral of motion. 
Indeed, to remain constant during the evolution of the system the Poisson bracket of 
a time-independent F with respect to the Hamiltonian H must vanish, 

{H,F} Po = X H (F) = 0. (2.10) 

The proof of this statement is straightforward upon considering the total time deriva- 
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2.1 Hamiltonian Mechanics 



tive of F = F(q\pi] t) 



dF 
~dt 



dF 
~dt 




By Hamilton's equations (2.6) and (2.5) this takes the form 



dF 
~dt 



dF 
~dt 



+ {H,F} Po , 



(2.11) 



where 



{H,F} Pe 



E 



dH dF dF dH 

dpi dq l dpi dq l ' 



is the Poisson bracket of H and F with respect to a canonical Poisson bivector P Q . 
In the case that F is time-independent, (2.10) follows immediately from (2.11). 

Hamiltonian systems play an essential role in classical and quantum mechanics and 
the theory of differential equations. They can be used to describe simple Newtonian 
systems such as the motion of a pendulum, or a heavy symmetrical top with a fixed 
lower point, and other dynamical systems such as planetary orbits arising in celestial 
mechanics. In this thesis, we will focus on a class of Hamiltonian systems whose 
Hamiltonian functions take the form of a natural Hamiltonian. To establish this 
distinction we will first note what it means in the case of the Lagrangian formalism, 
as this readily motivates the interpretation in the Hamiltonian formalism. 

In the context of the Lagrangian formalism, a natural Lagrangian is one which 
takes the form 



where gij denotes the covariant components of the associated metric tensor g, which 
is a function of the coordinates q l , and V(q l ) represents the potential energy, which 
characterises the interaction between the particles in a system. We remark that in 
physical terms the Lagrangian is defined as the difference between kinetic and poten- 
tial energy of a mechanical system, which explains the aforementioned terminology. 
With this identification, the first term of (2.12) can be interpreted as the kinetic en- 
ergy of the system. For the natural Lagrangian (2.12) the Euler-Lagrange equations 




(2.12) 
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2.1 Hamiltonian Mechanics 

(2.2) written in covariant form appear as 

H ■ Vj^'f H ; ' V r ( 2 - 13 ) 

where T l - k are the components of the Levi-Civita connection discussed in Section 
1.3. In the case that we have a vanishing potential, i.e. a system free of particle 
interactions, (2.13) takes the form of the geodesic equation, which is very familiar: 
the geodesic equation models the evolution of a system with n degrees of freedom 
and no external forces. Through the Legendre transformation (2.4) we then obtain 
the equivalent natural Hamiltonian, which assumes the form 

H(q\Pi) = b ij PiP3 + v(q% (2-14) 

which our discussion above makes it clear that the Hamiltonian can be interpreted 
as the total energy of our system, i.e. the sum of kinetic and potential energy, respec- 
tively. 

Neither the Lagrangian nor the Hamiltonian formulation present a simple system 
of equations to solve. Even though the Hamiltonian formulation reduces the Hamil- 
tonian system to a set of first-order equations, the fact remains that we still must 
solve a coupled set of non-linear ordinary differential equations. This prompts the 
need for methods by which we can complete the complicated task of integrating the 
equations of motion. In the case of systems with two degrees of freedom, we can 
focus on the method of separation of variables in the context of the Hamilton- Jacobi 
theory, which will allow us to link integrability of the equations of motion with the ex- 
istence of an associated valence two Killing tensor. This link is in part a result of the 
Arnold-Liouville theorem [4] which states that a Hamiltonian system with n degrees 
of freedom will be integrable by quadratures if it admits n functionally independent 
first integrals of motion in involution with the Hamiltonian. 
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2.2 Hamilton-Jacobi Theory 

The strategy behind the Hamilton-Jacobi theory is to rewrite Hamilton's equations, 
(2.5) and (2.6) in a revealing form through a choice of a new system of coordinates 
(see e.g. [39, 4]). In particular, we will require a canonical transformation of the 
coordinates, i.e. a transformation that leaves Hamilton's equations invariant. Indeed, 
a canonical transformation would allow us to obtain an equivalent formulation of 
classical mechanics just as in the previous section. Explicitly, if we change to a 
new system of coordinates (Q l , Pi) through a canonical transformation then the new 
Hamiltonian K must satisfy 

di _dK p= _8K 

To this end we can consider the action (2.1) as a function of the coordinates, i.e. 
independent of time. We will then study the change in the action S, defined in the 
previous section, from one path to a neighbouring path, which is readily calculated 

as 



8S 



dL 11 



oq 



U Ju \ d Q dtdq) q 



Since the path of the system satisfies the Euler-Lagrange equations (2.2), then 
the integral in bS vanishes. If we consider the canonical momenta coordinates (2.3) 
then we see that the above equation gives us, for an arbitrary number of degrees of 
freedom, 

< 2 - i5 » 

From the definition of the action (2.1) we also have the necessary link between this 
alternative formulation and the Lagrangian formulation. In particular we see that 
the time-derivative of the action S(q\t) is equal to the Lagrangian of Section 2.1. 
Namely, 

dS/dt = L. (2.16) 
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2.2 Hamilton-Jacobi Theory 



Combining this with (2.15) then yields the following relationship, 

dS{q\t) _dS dSdq i 
It ~~ ~dt + dq*H 

= 9t+p*?- ( 2 - 17 ) 

Now consider that we seek a canonical transformation. In this case by imposing 
that these new coordinates (Q l , Pi) also satisfy the principle of least action we will 
derive an analogous relationship which must hold with respect to the new Hamilto- 
nian. Thus we can take the difference between the new action and the original action, 
which we write as 

dF 

Pi q l - PiQ 1 = H - K + — , 

dt 

where F = F(t,q l ,pi,Q l , Pi) defines a generating function that characterises the 
canonical transformation used to define the new coordinates. We can determine the 
appropriate Legendre transformation in the above equation by first taking our new 
generating function to be G = F + Piq 1 , so that it depends on the old position 
coordinates and the new momenta coordinates. With respect to the canonically 
conjugate quantities q % and Pi, we then arrive at the following system of equations 

= v Qi = W K = H+ W (2 ' 18) 

The generating function will be chosen in such a way that the new Hamiltonian K is 
identically zero 1 (see e.g. [39]). In this case the above (final) equation (2.18) yields 
(replacing with G = S) 

83 hU°*)=0, (2,9) 



dt \ dq l 

which one can take as the Hamilton-Jacobi equation, though in this thesis, we will 
reserve this term for another version of the equation. 

We proceed a bit more in the derivation upon noting that the Hamilton equations 
(2.5) and (2.6) in the new coordinates become trivial, i.e. Q % = Pi = 0, so that 
Q % = d l and Pi — where d l and Cj represent an n-tuple of constants. Thus, we can 



lr This choice refers to a system which has thus been linearised. 
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2.2 Hamilton-Jacobi Theory 



write 



S = S(t,q l ;d), pi 



dS 
dq i 



dS 
dot 



Since the natural Hamiltonian (2.14) does not have any explicit time dependence, 
then we can let 



which makes it clear that the time dependence is trivial. In particular we can choose 
So(t; E) = —Et, which corresponds to the separation of the time variable since E 
only represents a constant. With this choice equation (2.19) then takes the form of 
what we will refer to as the Hamilton-Jacobi equation 



A complete integral of the HJ equation will then be of the form W = W(q l ; q) which 
satisfies the non-degeneracy condition 



By Jacobi's theorem once a complete integral is known then one can compute the 
motion of the system explicitly, i.e. determine the trajectories of the Hamiltonian 
system. And so finding a complete integral to the Hamilton-Jacobi equation affords 
us an alternative way to derive a solution to the associated Hamiltonian system. 

Theorem 2.22 (Jacobi). Let W (q l ; q) be a complete integral of the Hamilton-Jacobi 
equation (2.20) and let t and d 1 , . . . , d n ~ l be arbitrary constants. Then the functions 



S(t,q l ;ci) = S (t; q 



) + W{q\c i ), 




(2.20) 




(2.21) 



q l = q l (t] d,d l ) 



defined by the relations 



dW 



dW 



t-to 



BE 



dci 



i = 



l,...,n-l 
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2.2 Hamilton-Jacobi Theory 



together with the functions 



dW 



Pi = 



i = 



form a general solution of the canonical Hamilton equations, (2.5) and (2.6). 

This result must be considered in connection with the Arnold-Liouville theorem 
which states that a system in n degrees of freedom will be integrable by quadratures if 
it has n functionally independent first integrals in involution (see, e.g [4]). In our case 
of systems with two degrees of freedom, i.e. defined on the Euclidean plane, we have 
as a corollary to the Arnold-Liouville theorem that finding one first integral F inde- 
pendent of the Hamiltonian H renders the system integrable by quadratures. Other 
criterion and results which determine the integrability of a general n-dimensional 
pseudo-Riemannian manifold have also been established by Stackel in the 1890's [60], 
and Levi-Civita in the early 1900's [41]. 

A powerful way of integrating a Hamiltonian system is crafted through an ap- 
propriate canonical transformation that places the associated HJ equation into a 
separable form. In this case the HJ equation is said to be integrable via separation of 
variables with respect to the newfound separable coordinates, which we will denote 
as u l . In this situation the separation ansatz for integrating the HJ equation takes 
an additive form according to 



which is subject to the non-degeneracy condition (2.21). Any Hamiltonian system 
for which there exists such a system of separable coordinates u % on Ai that yields a 
complete integral of the form (2.23) is termed to be separable. Orthogonal separation 
of variables occurs in the case that the canonical transformation is a point transfor- 
mation that leaves the metric tensor g diagonalised with respect to the new separable 
coordinates. Such a Hamiltonian system is said to be orthogonally separable. In the 
next section we elaborate on the historical development of separation of variables on 




(2.23) 



E 2 . 
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2.3 Separability on the Euclidean Plane 

One of the primary results of relevance to the present study comes from Liouville in 
1846 [42]. Liouville studied Hamilton's equations for the motion of a particle on a 
curved surface which was under the influence of a time-independent potential. He 
found that if the metric and potential of a corresponding Hamiltonian function took 
a special (separable) form then the Hamiltonian system was solvable. 

Theorem 2.24 (Liouville). Let Ai be a two-dimensional manifold with local position 
coordinates (u,v) and corresponding canonical momenta coordinates (p u ,Pv)- If a 
Hamiltonian function is of the form 

H = (A(u) + B{v)Y l [| ( P l + pi) + C{u) + D(v)] , (2.25) 

where A(u),B(v),C(u) and D(v) are arbitrary smooth functions, then the Hamilto- 
nian system given by {Ai,X h, H} can be solved by quadratures. 1 

Through the equivalence of the Hamilton- Jacobi theory, Liouville's result also 
demonstrated that the associated HJ equation of (2.25) would be solvable under 
additive separation of variables. In covariant form, we can write the metric of the 
kinetic part of (2.25) as 

ds 2 = (A(u) + B(v))(du 2 + dv 2 ), (2.26) 

so that any potential which appears as in the Hamiltonian (2.25), or any metric of 
the form (2.26) is said to be in the Liouville form. 

In 1881, Morera proved the converse of Liouville's result [47]. This established 
that if a system defined by a natural Hamiltonian (2.14) was separable in the context 
of HJ theory then its metric g and potential V(q l ) would take the Liouville forms 
with respect to the separable coordinates. Nevertheless, Morea's equivalence did 
not provide a way of determining the separable coordinates, nor did it address the 
possibility that a Hamiltonian system could be separable in some other system of 
coordinates, with respect to which the metric and potential were not in the Liouville 

1 By quadratures, one means through algebraic means or by taking an integral. 
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form. Bertrand and Darboux, in 1857 [6] and 1901 [16], respectively, worked to 
address these concerns. 

Bertrand searched for natural Hamiltonian systems (2.14) that could admit a first 
integral of motion which took the form 

F{q\ Vi ) = fiT J '(g>^ + U(q l ), i,j = 1,2, (2.27) 

for some K and U in the coordinates q % . Such a first integral is at most quadratic 
in the momenta. Bertrand showed that subject to the forms (2.14) and (2.27), the 
vanishing of the Poisson bracket of H and F, namely {H, F}p = 0, imposed two 
conditions on the Hamiltonian system. The first of these is the Killing tensor equation 

[g,K] = 0, (2.28) 

where [ , ] denotes the Schouten bracket presented in Section 1.3. The second is the 
compatibility condition 

d (KdV^j = 0, (2.29) 

which in (local) coordinates yields what is called the Bertrand- Darboux PDE (see 
e.g. [55]). The compatibility condition appears as the integrability condition imposed 
on U(q l ) in (2.27) by requiring that H and F be in involution, namely 

dU = KdV, (2.30) 

where K is the (l,l)-tensor defined by K := Kg 1 , i.e. K % j := K %i ggj. Together, 
the Killing tensor equation (2.28) and (2.30) are equivalent to the vanishing of the 
Poisson bracket of H and F subject to (2.14) and (2.27). 

In 1901, Darboux solved the linear, second-order PDE admitted by (2.29) using 
the method of characteristics [16], by first simplifying the PDE through a coordinate 
rotation and translation. This essentially had the effect of transforming the associated 
Killing tensor in (2.27) to a canonical form, which he could then diagonalise to 
derive the appropriate separable coordinates (u, v) [55]. In a systematic way, Darboux 
arrived at a potential satisfying (2.29) with respect to a canonical Killing tensor whose 
normal eigenvectors generated elliptic-hyperbolic coordinates. 
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Indeed, the problem that Darboux solved was to determine the most general 
potential which would admit orthogonal separation of variables with respect to an 
elliptic-hyperbolic coordinate system. This method is very familiar in classifying 
systems of superintegrable potentials, where one requires that a system be separable 
with respect to some particular orthogonal coordinate systems and then solves for the 
most general potential compatible via (2.29) (see Section 3.1). The complicated part 
of the above procedure rests in solving (2.29), which yields a non-linear second-order 
PDE when expressed in local coordinates. 

In 1934, a new approach to this problem was presented in a famous paper by 
Eisenhart [20] who provided an intrinsic characterisation of orthogonally separable 
Hamiltonian systems. Eisenhart worked to extend the results of Stackel [59], who 
studied the orthogonal separability of the HJ equation. Stackel established that a 
necessary condition for a system defined by a natural Hamiltonian to be orthogonally 
separable was that it admits n — 1 quadratic first integrals of the form (2.27), all 
functionally independent and in involution with H. Eisenhart observed that a smooth 
function F G T*(M) in the restricted form 

F{q\ Vi ) = \K l \q i ) VlV j i,j = 1,2, (2.31) 

would be a first integral of the geodesic Hamiltonian, 

Htf,Pi) = ±g ij p i p j , (2.32) 

if and only if the functions K % i in (2.31) were the components of a characteristic 
Killing tensor field K e K?{M). This result is readily established upon computation 
of the Poisson bracket of F and H in the assumed forms (2.31) and (2.32), respectively. 
The pivotal result Eisenhart established in the theory of orthogonal separation of 
variables is given by 

Theorem 2.33 (Eisenhart). The Hamiltonian system defined by a geodesic Hamil- 
tonian (2.32) is orthogonally separable if and only if it admits n — 1 functionally 
independent first integrals of motion taking the form (2.31), such that 

(1) all corresponding Killing two-tensors have real and pointwise distinct (almost 
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everywhere ) eigenvalues, 

(2) all corresponding eigenvector fields of the Killing two-tensors are normal, 

(3) the Killing two-tensors defined by the n — 1 first integrals all have the same 
eigenvectors. 

A Killing tensor satisfying (1) and (2) is called a characteristic Killing tensor. 
Eisenhart's results demonstrated that Killing tensor fields would play a crucial role 
in determining the orthogonal separability of a Hamiltonian system. In 1997, Benenti 
[5] extended this result to Hamiltonian systems defined by a natural Hamiltonian, 
i.e. with non-vanishing potential. 

Theorem 2.34 (Benenti). The Hamiltonian system defined by a natural Hamiltonian 
(2.14) is orthogonally separable if and only if there exists a characteristic Killing two- 
tensor K such that 

d (Kdv) = 0. 

The requirement of n — 1 first integrals in Eisenhart's theorem is reduced to the 
existence of a single characteristic Killing tensor in Benenti's theorem. To determine 
such a Killing tensor, one can start with the n — 1 Killing tensors afforded by the 
n — 1 first integrals of motion, and take a linear combination of these Killing tensors. 
If we let Ki, ... , K n _i be as in Theorem 2.33, then adding the metric tensor g of Ai 
will yield a basis generating an n-dimensional vector subspace of K?(M). A general 
Killing tensor in this space 

n-1 

K = g + J2 K - 

i=i 

is then guaranteed to still have real and distinct eigenvalues and the same eigenvectors 
admitted by the individual K j, i — 1, . . . ,n — 1. The condition that the eigenvectors 
be normal, that is, 

Ei A dEi = 0, % = 1, . . . , n (no summation), (2.35) 
where Ei denote the eigenforms of K, then implies that the eigenforms generate n 
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foliations 1 . These will consist of (n — l)-dimensional hypersurfaces orthogonal to 
the eigenvectors of each Killing tensor (see [8, 9] for details). This is the geometric 
construction of the orthogonal coordinate web. Such a web defines the separable 
coordinates with respect to which an associated HJ equation separates. 

These two theorems establish the explicit role characteristic Killing tensors play 
with regards to the HJ theory of orthogonal separation of variables, and essentially 
provides a link between their algebraic and geometric properties. In the next section 
we will explore this link through the theoretical framework of the invariant theory of 
Killing tensors with an emphasis on defining joint invariants, as first introduced in 
[56], and later utilised in the theory of superintegrable Hamiltonian systems in [3]. 

2.4 Invariant Theory of Killing Tensors 

We begin this section with a brief historical account on the emergence of Killing 
tensor fields in mathematics. For a more in-depth account of this history, we direct 
the reader to [27], on which this introduction is based. 

In the 1880s, Wilhelm Killing began an extensive study of non-Euclidean geometry 
in which he sought to develop a theory of space forms, which in modern language is a 
complete Riemannian manifold A4 with constant curvature. Killing began his study 
analytically by studying the behaviour of infinitesimal motions of an n-dimensional 
continuous manifold of points (xi, . . . , x n ) with n degrees of freedom. In his attempt 
to deal with all possible space forms, Killing imposed an (unconventional) condition 
on his infinitesimal motions which lead to the implication that they formed a finite- 
dimensional Lie algebra, a theory mostly unknown to him at the time. This naturally 
brought Killing into contact with Sophus Lie, and through his and Lie's results on 
transformation groups he succeeded in proving that a proper space form would have 
degree n(n+l)/2 and admit a Riemannian metric. From the fact that an infinitesimal 
motion ought to leave the metric invariant, Killing then arrived at the equations 

1 In the case of E 2 , we have that (2.35) is always satisfied, by dimensional considerations, and 
so that the eigenvectors of Killing tensors defined on E 2 are always normal. 
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defining what we now call a Killing vector field 



C K g = o, 



(2.36) 



where Ck denotes the well-studied Lie derivative along the vector field K, and g is 
the metric tensor of the Riemannian manifold M.. Indeed, we can use this Killing 
vector field to define a function on the cotangent bundle T*(Ai), namely K l pi, and 
(2.36) then implies that this function will be a first integral of the geodesic flow on 
T*(Ai) with respect to the Riemannian metric g. 

2.4.1 Preliminaries 

The Killing fields with which Killing himself worked were Killing vector fields, which 
act as the infinitesimal generators of isometries on a Riemannian manifold. This 
is made obvious by considering equation (2.36), where we observed that the vector 
fields K admitted by this equation were those which preserved the metric g. Killing 
tensor fields of higher valence, p > 1, provide information regarding quadratic, cubic, 
and higher-order first integrals of the Hamiltonian geodesic flow. With this in mind 
we begin our study of the vector space in (Ai,g) formed by Killing tensors of the 
same valence p, which we will denote as KP{M). The fact that this collection of 
tensors forms a vector space follows from the bilinear nature of the Schouten bracket 
(see Section 1.3). 

If Ai is a space of constant curvature then the dimension of 1C p (M. n ) will be 
maximal, and given by the Delong-Takeuchi-Thompson (DTT) formula [18, 61, 64] 



In this situation, we see that the vector space 1C p (Ai) is determined with respect 
to d arbitrary parameters. In other words, an element of IC p (Ai), namely a Killing 
tensor with fixed valence p, can be viewed as being an algebraic object in a vector 
space. This identification provides the pivotal link between vector spaces of Killing 
tensors defined on pseudo-Riemannian spaces of constant curvature and the classical 
invariant theory of vector spaces of homogeneous polynomials (for more details on 




(2.37) 
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the latter see e.g. [49]). 

The focus of this thesis is the vector space of Killing two-tensors defined on the 
Euclidean plane, namely /C 2 (E 2 ) and products of this space, e.g. /C 2 (E 2 ) x /C 2 (E 2 ). 
The DTT formula yields that the dimension of /C 2 (E 2 ) is six, since n = p = 2. An 
alternate way of arriving at this number is to formally derive the solution to the 
Killing tensor equation (2.28) in the Euclidean plane, which we do in Appendix B, 
and count the total number of integration constants in the end. Additionally, this 
approach yields that the general solution to (2.28) in E 2 , with respect to Cartesian 
coordinates, is a Killing two-tensor K G /C 2 (E 2 ), with components given by 

K ll = fa + 2fay + fay 2 , 

K 12 - fax - fay - faxy, (2.38) 
K 22 = fa + 2 fax + fax 2 , 

where fa,i = 1,. . . ,6 are arbitrary parameters which appear as constants of inte- 
gration in solving (2.28). And so indeed, we can identify the vector space of Killing 
two-tensors defined on E 2 with IR 6 . Another result which will be frequently made 
use of throughout this thesis is the following lemma, independently arrived at by 
[18, 61, 64] in the early 1980's: 

Lemma 2.39. Any Killing tensor of valence p defined on a constant curvature 
pseudo-Riemannian manifold (Ai,g) can be expressed as a sum of symmetrised tensor 
products of a basis of Killing vectors on Ai . 



as 



We make use of the above lemma to write our general Killing tensor K e /C 2 (E 2 ) 
K = A ij Xi Xj + B'Xi R + CRQ R, (2.40) 



where Xi denote the translational basis Killing vectors and -R denotes the rotational 
basis vector. In this case the A % \B % and C define the Killing tensor parameters, 
which we can identify with (2.38) by letting 



A-= ? " , 5*= Z , C = fa, (2.41) 
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which will be subject to the symmetry property 

A ij = A (ii)_ 

Furthermore, we will take as a basis for the space of Killing vectors on E 2 the usual 
vectors given in Cartesian coordinates, {q l ,q 2 ) = {x,y), namely 

d d d d 

X 1 = ^-, X 2 = ^-, R =y lL- x ° (2.42) 
ox oy ox oy 

The components of the general Killing tensor (2.40) with respect to the basis (2.42) 
are then expressible in a compact way as 

K ij = A ij + 2e {i e B j) x e + Ce i m e j k x m x k . (2.43) 

One can check (see Section 1.3) that indeed, subject to (2.41) this reproduces the 
general Killing tensor components defined in (2.38). In this compact presentation, it is 
straightforward to compute the transformation law for the parameters fy, i — 1, . . . , 6 
under the action induced by the isometry group SE(2) O E 2 which acts transitively 
on the Euclidean plane via 

q { + 5\ (2.44) 

where A l j G 50(2), 8 G IR 2 , so 5 l = Pi,i = 1,2. This computation is carried out 
using the compact notation in Appendix A.l. For clarity we will now provide the 
explicit computation of the induced action. 

We seek to derive the transformation law for the parameters in /C 2 (E 2 ) under the 
action of SE(2), which we can obtain by making use of the usual transformation rules 
for tensor components (mentioned in Section 1.3). In particular, a Killing two-tensor 
transforms according to the tensor transformation rules via 

K ij = A i iA j kK tk , (2.45) 

where the A*j represent the Jacobian matrices (1.6) defined in Section 1.3. The action 
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of SE(2) E 2 is identified by 

q 1 \ ( cos 9 — sin 9 
q 2 J \ sin (9 cos 9 

so that identifying the components of the Jacobian via (1.6), namely 




(2.46) 



A' t = g, (2.47) 



we see that indeed this is given by a rotation matrix: 
A = 



dq l /dq x dq^/dq 2 \ ( cos 9 —sin 9 
dq 2 /dq l dq 2 /dq 2 ) \ sm.9 cos 9 



Thus, computationally, equation (2.45) amounts to (appropriately) multiplying the 
2 x 2-matrix defining the Killing two tensor K % i with two rotation matrices: 

K = \ T K\, 

where we see that we must multiply on the left by the transpose of the Jacobian 
matrix in order to carry out the proper tensor index transformation. Carrying out 
this calculation, we obtain in a simple manner that the components of K % i transform 
according to 

K u = -K u cos 2 9 + 2K 12 cos 9 sin 9 + K 22 sin 2 9, 

K 12 = -K 11 cos 9 sin 9 + K 12 (cos 2 9 - sin 2 9) + K 22 cos 9 sin 9, (2.48) 
K 22 = K 11 sin 2 9 - 2K 12 cos 9sin9 + K 22 cos 2 9. 



Using the compact formulae (2.43), the non-transitive action of SE(2) /C 2 (E 2 ) 
is elegantly derived by substituting into (2.40) with the basis vector transformation 
rules (A.l), which yields equation (A.3) 1 . Since we work in E 2 we can get by without 
the compact notation, and so we continue with the above procedure for deriving the 

1 See Appendix A.l for this computation. 
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invariants as follows. First, we must substitute with the polynomial components of 
K l 3 given by (2.38) into the transformed components (2.48) and then set this equal 
to K l i written in transformed parameters, i.e. we arrive at the following system of 
equations (we have let (q 1 ^ 2 ) = (x,y) for clarity) 

fo + my + Pstf = + Wav + fay 2 ) cos 2 e 

+ 2(fo — fox — foy — foxy) cos 9 sin 9 
+ (fo + 2f3 5 x + fox 2 ) sin 2 9 
fo - fox - foy - foxy = -(/3i + 2foy + foy 2 ) cos 6 sin 9 

+ (fo ~ fox - foy - foxy)(cos 2 - sin 2 0) 
+ (fo + 2 fox + fox 2 ) cos 9 sin 9 
fo + 2 fox + fox 2 = (fo + 2 foy + foy 2 ) sin 2 6 

— 2(fo — fox — foy — foxy) cos 6 sin 6 
+ (fo + 2 fox + fox 2 ) cos 2 9 

Next, we must substitute for x and y according to (2.46), after which we can equate 
coefficients of like terms. This then yields the following explicit transformation laws 
for the parameters of /C 2 (E 2 ): 

fo = focos 2 6 - 2/3 3 cos 6> sin 6> + fosm 2 6 - 2p 2 focos6, 

-2p 2 fo> sin 9 + fopj, 
fo = fosm 2 6 - 2/3 3 cos6>sin6>, +focos 2 6 - 2p 1 focos6, 

+2 Pl fosm6 + fop 2 , 
fo = (ft-^^sin^cose + ^^^-sin 2 ^, (2.49) 

+ (pifo + P2fo) cos 9 + (pifo -p 2 fo) sintf - fopip 2 , 
fo = focos9 + fosin9 - fop 2l 
fo = fo cos 9 - fo sin 9 - fop ly 

fo = fo- 

This brings to light our first (algebraic) invariant, namely Ai = fo. 1 In order to 
1 See Appendix A A. 2 for the compact form of these six parameters under the group action. 
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derive the remaining fundamental invariants we use the method of moving frames. 
2.4.2 The Moving Frames Method 

In what follows we provide a synopsis of the moving frames method and its role in 
the invariant theory of Killing tensors. For a more complete account we direct the 
reader to the wealth of literature that has nourished the development of this theory, 
namely [45, 17, 1, 3, 13, 12, 15, 29, 30, 31, 32, 55, 56, 70, 73] and references therein. 
In particular, [45] pioneered the study of invariants under the isometry group acting 
on K p {Ai) into the theory of Killing tensors as a way of classifying the orthogonal 
coordinate webs in the Euclidean plane. 

As noted in Section 2.3, it has been shown (see e.g, [20, 5, 31]) that an element 
K e K?(M) with real and distinct eigenvalues and normal eigenvectors generates 
an orthogonal coordinate web with n = dim(A^) foliations whose leaves are n — 1 
dimensional hypersurfaces orthogonal to the eigenvectors of K. In this situation 
Cartan's geometry [11] can be combined with the study of principal fiber bundles to 
provide a framework for the invariant theory. For a detailed account of the underlying 
ideas of the generalised moving frames method we direct the reader to [49], from 
which this brief overview is based on. In this section, we focus on demonstrating 
the recursive approach to constructing the moving frame map, developed by Kogan 
in 2003 [38]. We begin with some preliminary definitions mostly compatible with 
[49, 56, 55, 45]. 

Let G be an r-dimensional (Lie transformation) group acting smoothly on an 
n-dimensional manifold Ai with s-dimensional orbits. Representative points in each 
group orbit can be chosen to depend continuously on the orbits provided that the 
group acts regularly 1 . A (local) cross-section is an (n — s)— dimensional submanifold 
K C Ai that intersects each group orbit transversally 2 and at most once. In general, 
a local cross-section passing through any point in p 6 Ai can be constructed provided 

1 Regularity means that all the orbits of the group action on Ai have the same dimension and 
each point p £ M has a system of arbitrarily small neighborhoods whose intersection with each 
orbit is a pathwise connected subset of the orbit [49] 

2 Two submanifolds K,N are said to intersect transversally at a point p G K D N if T p (K) D 
T p (N) = {0},i.e. if they share no common non-zero tangent vectors. 
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that the the group acts regularly on M.. 

Choosing a cross-section amounts to fixing a moving frame, i.e. identifying a map 
/ : K, 2 (Ai)/G — > G, where we focus on the quotient space, considering G as the 
subspace G = {cg\c e R} generated by the metric, which in our case is trivial since 
it provides no information regarding separable coordinate systems [45]. A moving 
frame is formally defined as a smooth, G-equivariant map 1 p : M. — > G whose exis- 
tence is guaranteed in a neighborhood of a point p G M. provided that G acts freely 2 
and regularly near p [49]. To construct a moving frame, one utilises Cartan's nor- 
malisation method, whereby the elements of the cross-section can be interpreted as 
canonical forms for general elements in the underlying manifold Ai. The normalisa- 
tion procedure amounts to choosing local coordinates g — (gi, . . . , g r ) on G near the 
identity element, identifying explicit formulae for the group transformations in the 
coordinates, and then equating the first r components of the formulae in the previous 
step to given constants. The Implicit Function theorem then guarantees that such a 
system of equations is locally soluble. 

We shall illustrate the method explicitly by taking G to be the isometry group of 
the Euclidean plane, which acts as an automorphism, i.e. Killing tensors are mapped 
to Killing tensors in /C 2 (E 2 ), and so preserves the geometry of the vector space. 
Explicitly, we take G = SE(2), the special Euclidean group, which acts freely and 
regularly in the vector space of Killing tensors [56] . 

Recall that SE(2) consists of the rigid body motions, translations and rotations, 
and so it can be identified as the semi-direct product of the special orthogonal group 
50(2) with the group of translations on E 2 . The recursive version of the moving 
frames method to derive invariants then proceeds in two steps. First, we compute 
a set of fundamental invariants under the subgroup of translations and choose an 
appropriate cross-section to obtain a set of translational invariants. Second, we use 
these translational invariants as new coordinates in the vector space and compute 
the action under the subgroup of rotations, 50(2). After choosing an additional 

1 A map 7 : A4 —> G is G-equivariant if for g G G, p € A4, ~f(g ■ p) — j(p) ■ g^ 1 with respect to 
the actions of G on M. and on itself by right multiplication. 

2 A transformation group G acts freely provided that the isotropy group, G p = {g G G\g ■ p = p} 
where p 6 A4 is trivial: G p = e, for all p € Ai, where e denotes the identity element of G. In other 
words, only the identity clement fixes any p € M. 
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cross-section, or normalisation equation, we then arrive at a final set of invariants 
under the action of the full isometry group [38]. 

The action of the subgroup of translations amounts to substituting with A\ = 5\ 
in (2.44) and (2.49), where 5^ denotes the Kronecker delta. In other words, the 
group of translations acts on E 2 via q l — > q l + 5 l and so we find that this induces the 
transformation 1 , 

01 = 0i + 20 4 p 2 + 6 pl 

02 = 02 + 205P1+M, 

03 = 03 - 04P1 - 05P2 ~ PePl&i, ^ 
04, = 04 + 06P2, 

05 = 05 + 06Pl, 

06 = 06- 

We then choose our cross-section for the restricted group action via 

04 = 5 = 0. (2.51) 



The moving frame map for the restricted group action is then defined by the following 
normalisation equations 

Pi = -04/06, P2 = -05/06 (2.52) 

We then substitute these relationships into the remaining four equations (2.50) and 
thus arrive at the following fundamental (translational) invariants of /C 2 (E 2 ) 

h = 0106 ~ 01 

h = 0206 ~ 01 (2.53) 

h = 0306 + 0405, 
h = 06- 



Next, we determine the action of the subgroup of rotations SO (2) on the space of 
the four invariants (2.53). Since 5*0(2) acts on E 2 via x — > A l jx\ then we can use 

1 See Appendix A. 2, eq. A. 3 for the compact form of this transformation. 
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the same technique as above to determine the transformation formulas induced by 
this action. Thus, we find that 5*0(2) acts on /C 2 (E 2 ) according to 1 



= 1 cos 2 9 + 2 sin 2 9 - 20 3 sin 9 cos 9, 

= 0i sin 2 9 + 2 cos 2 9 - 20 3 sin 9 cos 9, 

03 = {0i — 02) sin 9 cos 9 + /^(cos 2 9 — sin 

04 = 0a cos 9 + 5 sin 9, 

R. — R. uir, ft — R. nr^a ft 



(2.54) 



05 = 04 sin 9 — 5 cos 9, 

06 = 06- 



Therefore, taking into consideration (2.53) and (2.54), we find that 

I 1 = - (cos 90 A + sin 60 5 ) 2 + (cos 2 90 1 + sin 9 (sin 90 2 - 2 cos 60 3 )) 6 , 
J 2 = - (sin 904 - cos 90s) 2 + (sin 2 90 1 + cos 9 (cos #/3 2 - 2 sin 0/3 3 )) /3 6 , 
7 3 = \ sin 20 (-/3 2 + 01 + (& - 2 ) 6 ) - cos 20 (0 4 5 + 3 6 ) , 
J 4 = h- (2.55) 

We can then obtain the final normalisation equation by taking .Z3 = 0, which yields 

n 1 ( 2(0 4 5 + 3 6 ) 

arctan 



2 y Wl -p 2 )-pl + pl 

Substituting this into (2.55), we recover the SE{2) invariants 

Ai = 6 , 

A 2 = 6 (0 1 + 2 )- 0t-0l (2.56) 
A 3 = (06(0i ~ 02) ~ 01 + 0lf + 4 (060s + 040 5 ) 2 , 

established first in [70] and later in [45, 56, 55]. As was shown first by Winternitz 
and Fris in 1965 [70] and again in 2002 by McLenaghan et al. [45] via the invariant 
theory of Killing tensors, the invariants A x and A 3 can be used to classify the orbits 



1 See Appendix A. 2 eq. A. 4 for compact formulae. 
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of JC 2 (E 2 )/SE(2) as follows 



Elliptic-hyperbolic : Ai 7^ 



A 3 ^0, 
A 3 = 0, 



Parabolic : A 



1 — 







Polar : Ai ^ 



Cartesian : A 



1 — 







A 3 = 0. 



The final tool we require from the invariant theory comes from considering the 
group action on the product space obtained by taking two copies of the vector space 
of Killing tensors: /C 2 (E 2 ) x /C 2 (E 2 ). This action leads to the notion of joint invariants 
of Killing tensors, the topic of discussion featured in our next section. 

2.4.3 Joint Invariants of Killing Tensors 

We now present the geometric properties of joint invariants of Killing two-tensors (to 
be defined below) as first introduced by Smirnov and Yue in 2004 [56]. Our notation 
will be compatible with Adlam et al. [3], who extended the study established in 
[56] and demonstrated its applicability to superintegrable systems. We specialise this 
development for the isometry group which acts on the product space /C 2 (E 2 ) x/C 2 (E 2 ), 
where we begin with the joint invariants of non-degenerate orbits of the orbit space 



A three-dimensional orbit is defined to be non- degenerate provided that, along the 
orbit, the invariant given by k 2 = y/A^/(A[) 2 , is non-vanishing. The action SE(2) O 
/C 2 (E 2 ) x /C 2 (E 2 ) for which k 2 7^ is likewise termed non-degenerate. The geometric 
interpretation of the invariant k 2 is the (half) distance between the foci of an elliptic- 
hyperbolic coordinate web generated by the normal eigenvectors of a Killing two- 
tensor, and so it follows that such Killing tensors will have a correspondence with 
the non-degenerate orbits [45, 55, 70]. 

In this situation, the isometry group SE(2) acts on each copy of /C 2 (E 2 ) with 
three-dimensional non-degenerate orbits. If we denote the parameters of each vector 
space by ai,(3i,i = 1,...,6, respectively, then the group action is induced by the 
corresponding transformation laws given by (2.49). Explicitly, we have that the 



2 (E 2 ) x /C 2 (E 2 ))/S£(2) studied in [3]. 
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parameters ojj, i — 1, . . . , 6 in the first copy of /C 2 (E 2 ) transform as 

5ii = oil cos 2 9 — 2a 3 cos 9 sin 9 + a 2 sin 2 — 2p 2 a 4 cos 

-2p 2 «5 sin6> + "6P2> 
a 2 = on sin 2 9 — 2a 3 cos 9 sin# + a 2 cos 2 9 — 2pia$ cos 

+2pia4 sin + oiqPi, 
&3 — («i — « 2 ) sin0 cos0 + a 3 (cos 2 — sin 2 0) (2-57) 

+ (pi« 4 + P2«5) cos + (pia 5 - p 2 a 4 ) sin - a 6 pip 2 , 
a 4 = OJ4COS0 + oj5sin0 — a 6 p 2 , 
($5 = a 5 cos# — a 4 sin# — a 6 £> 1? 

"6 = "6, 

and the equivalent relationship in terms of the parameters (3i,i = 1, ... ,6, for the 
second copy, 

f3 l = /3 1 cos 2 0-2/3 3 cos0sin0 + /3 2 sin 2 0-2p 2 /3 4 cos0 

-2p 2 f3 5 sin# + /3 6 p 2 , 
/3 2 = / 0isin 2 0-2/9 3 cos0sin0 + /9 2 cos 2 0-2pi/9 5 cos0 

+2pi/3 4 sin^ + /3 6 p 2 , 
/3 3 = (/3 1 -/3 2 )sin0cos0 + /3 3 (cos 2 0-sin 2 0) (2.58) 

+(pi& + P2P5) cos 9 + (pi/3 5 - p 2 f3 4 ) sin 9 - /3 6 PiP2, 
(3 4 = /3 4 cos0 + /3 5 sin0 - (3 6 p 2 , 
h = fecosO - fasinO - fop!, 
h = As- 

Additionally the conditions 

2 _ a/ (aj - ol\ + a 6 (a 2 - cti)) 2 + 4(q 6 a 3 + c^c^) 2 
fc i — r 1 U, 

,2 _ v/(/3 2 - /35 2 + ggg - ft)) 2 + 4(/3 6 /? 3 + ftft) 2 , n 

K 2 — o T U 

hold true. 

If we denote the six dimensional space defined by the parameters ojj, i = 1, . . . , 6 by 
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A ~ M 6 , then we can define a, joint invariant [56] of the product space /C 2 (E 2 ) x/C 2 (E 2 ) 
ClS db function J:/1xj4->1 that satisfies 



under the transformation laws induced by the isometry group SE(2). Six joint in- 
variants of the action SE(2) /C 2 (E 2 ) x /C 2 (E 2 ) then immediately follow from (2.56), 
namely 



a; = a 6 , 

A' 2 = a 6 (ai + Oi 2 ) - a\ - a 2 , 

A' 3 = (ae(ai - a 2 ) - af + a 2 ) 2 + 4(a 6 a 3 + a 4 a 5 ) 2 , 

A' 4 = 

A' 5 = /3 6 (/3 1 + /3 2 )-/3 4 2 -/3 5 2 , 

A' 6 = ( / 3 6 ( / 3 1 - / 3 2 )- / 3| + / 3|)2 + 4(/3 6 /3 3 + /3 4 /3 5 ) 2 . 



We remark that since the joint invariants are obtained via the moving frames method 
then they are functionally independent. Further, any analytic function of the in- 
variants (2.61) will also be a joint invariant of the non-degenerate group action 
SE{2) /C 2 (E 2 ) x /C 2 (E 2 ) [49]. The Fundamental Theorem on invariants of regular 
Lie group actions tells us that we should obtain a total of 9 fundamental invariants, 
which we arrive at by subtracting the dimension of the orbits, 3, from the dimension 
of the product space /C 2 (E 2 ) x /C 2 (E 2 ), 12. Thus, we are in search of three more 
fundamental invariants. 

These remaining joint invariants can be derived by employing geometric consid- 
eration. Recall that each element of a non- degenerate orbit of the action /C 2 (E 2 ) x 
/C 2 (E 2 ) corresponds to a Killing tensor whose normal eigenvectors generate an elliptic- 
hyperbolic web. The foci represent the singular points of this coordinate web, which 
are characterised, following [5], as the points where the eigenvalues of the associated 
Killing tensor are equal. Establishing this degeneracy, one can solve for the explicit 
form of the singular points in terms of the arbitrary parameters, from which we have 
that the coordinates for the two foci, denoted S\, S2, in the elliptic-hyperbolic web 



J = F(a u . . . ,a 6 ,/3i, . . . , /3 6 ) 
= F(di, . . . ,d 6 ,/5i, • • • , As), 



(2.59) 
(2.60) 
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are (see [45]) 

(xi,Vi)si = 




(2.62) 

(X2, V2)s 2 = 

/ / — \ V 2 / / — 

-Ai 1 ( V A 6- CT 1 \ -§A 1 ( V A 6+^l 

ft ft I 2 j > ft ft I 2 

where <7i = /3| — /?f + /^(/^ — fa) and A' 6 is as denned in (2.61). Considering the 
product space /C 2 (E 2 ) x /C 2 (E 2 ), we have another set of foci, (x^, 2/3)53 and (#4, 2/4)54, 
appropriately defined with respect to the parameters cti,i = 1,...,6, o\ in terms 
of «j, and A3. In [3], this identification motivates the authors to take the action 
SE{2) O /C 2 (E 2 ) x /C 2 (E 2 ) as the free and regular action SE{2) E 2 x E 2 x E 2 x E 2 , 
and hence, use the Weyl theorem on joint invariants to conclude that the square of 
the distances between the foci will be a joint invariant of the preceding action. Thus, 
in the case of a non-degenerate action, we can take the remaining three invariants as 

A' 7 = d 2 (S 2 ,S 3 ) = (x 2 -x 3 ) 2 + (y 2 -y 3 ) 2 , 

A' B = d 2 (S 1 ,S 3 ) = (x 1 -x 3 ) 2 + (y 1 -y 3 ) 2 , (2.63) 

A' 9 = d 2 (S 2 , S 4 ) = (x 2 - x 4 ) 2 + (2/2 - 2/ 4 ) 2 , 

where (x^, yi),i — 1, . . . , 4 are given by (2.62) in terms of fa and analogously in terms 
of aii. Note that we require the quadrilateral S1S3S4S2 to be rigid so that we need to 
only specify the distance of one of the diagonals and the other can be determined as a 
function of these nine invariants. Other geometrically motivated joint invariants can 
be established in this manner as well, which would be obtained as functions of the 
invariants A[ — A' 9 , such as the areas within the quadrilateral whose vertices are given 
by the four singular points — see Figure 2.1. In the next chapter, we characterise the 
case when one of the Killing tensors K £ /C 2 (E 2 ) belongs to a degenerate orbit of 
the orbit space IC 2 (E 2 )/SE(2). 
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Chapter 3 

Superintegrable Systems on the 
Euclidean Plane 



We recall the significance of the first integrals of motion for a Hamiltonian system. 
These are smooth functions in the generalised canonical coordinates that remain 
constant along the Hamiltonian flow, or orbits, of a system. As such, they can be 
identified as constants of motion for a system which elicits many familiar examples of 
first integrals. Indeed, the total energy, angular momentum, and linear momentum 
are all examples of constants of motion that provide physical insight about a Hamil- 
tonian system: in particular, each of these is an example of a conserved quantity 
that corresponds to a symmetry of the Lagrangian. A systematic way of deriving 
these quantities from symmetry considerations is the heart of Noether's theorem (see 
e.g. [23]) where, for instance, conservation of energy is an artifact of the Lagrangian 
being invariant under changes in time. In brief, we see that the first integrals of 
motion help to establish physical properties of a Hamiltonian system by signifying 
the existence of structure and symmetry. 



3.1 Development of Superintegrable Systems 

A Hamiltonian system with n degrees of freedom is said to be completely integrable 
if it admits n globally defined, functionally independent first integrals of motion in 
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involution with respect to the Poisson bracket (see Section 1). In the case that 
such a system admits n + 1 such first integrals of motion the system is said to be 
minimally superintegrable, whereas it is termed to be maximally superintegrable if it 
admits 2n — 1 such first integrals. In some cases we find that a system is solvable 
via orthogonal separation of variables in more than one system of coordinates, which 
qualifies such systems as being multi- separable. For the purpose of this study, we 
focus on superintegrable systems that exhibit this property of being multi-separable. 

As discussed in Section 1.2, we can study the first integrals of natural Hamiltonian 
systems 

H(q i ,p l ) = lg i ' PlPj + V(q i ), (3.1) 

defined on flat pseudo-Riemannian manifolds which take a form quadratic in the 
momenta according to 

F(q\ Pi ) = \K^^)p lVd + U(q l ), i,j = 1,2, (3.2) 

which will afford orthogonal separation of variables in the associated Hamilton- Jacobi 
equation in the case that K l i is a characteristic Killing two-tensor satisfying the 
compatibility condition 

d (KdV^j = 0. (3.3) 

This is the content of Benenti's theorem (see Section 2.3). 

We can generalise this development to the study of first integrals polynomial in 
the momenta of degree p upon considering Killing tensors of valence p (for a recent 
account, see e.g. [28] ) 1 . Indeed, the search for first integrals of motion in higher 
degrees of the momenta has been extensively studied for over a century, with Drach 
initiating the study of first integrals cubic in the momenta in 1908 [19, 52]. In 2002, 
Gravel and Winternitz considered superintegrable systems on E 2 characterised by 
one linear and one (non-trivial) cubic first integral [25]. This work was continued 
by Gravel who then characterised systems with one quadratic first integral (derived 
from the system admitting separability in Cartesian coordinates), and one cubic first 

1 Though, in studying Killing tensors of valence p > 2, we lose the interpretation of a system 
being orthogonally separable. 
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integral [24]. The origins of the invariant theory approach can be traced back to 
Winternitz and Fris in 1965 [70], in which the authors computed the invariants of 
second-order symmetries of the Laplace equation defined on E 2 under the action of the 
isometry group, and then used these invariants to classify the orthogonal separable 
webs on E 2 . In 2002, this geometric approach resurfaced with McLenaghan et al. in 
[44], who presented it in the language of the invariant theory of Killing tensors. 
In 2008, Adlam et al. [3] extended the theory constructed in [45] and derived an 
invariant characterisation of the superintegrable potential admitted by the Kepler 
problem through the use of invariants and joint invariants of Killing two tensors that 
determined first integrals according to (3.2). 

In this thesis, we continue the work of classifying superintegrable potentials on 
the Euclidean plane by providing a new perspective on the ideas developed in [3]. 
In particular, we investigate the joint invariants of a degenerate and non-degenerate 
orbit of the orbit space (/C 2 (E 2 ) x K? {¥?))/ SE{2) admitted by a Killing tensor whose 
eigenvalues generate an elliptic-hyperbolic web in the canonical position, and a sec- 
ond Killing tensor whose eigenvalues generate a polar web in the non-canonical po- 
sition. To this end we begin by providing a joint invariant characterisation of the 
Smorodinsky- Winternitz system, 

H = Upl+pl)-" 2 (* 2 + y 2 ) + ^ + ^ (3-4) 

x y 

where a, (3, and u are arbitrary parameters. Our goal is to use joint invariants 
of Killing two tensors K defined on E 2 [56] to provide a link between the arbitrary 
parameters of the SW potential and the parameters which characterise the associated 
vector space of Killing tensors /C 2 (E 2 ). 

Further, this will provide motivation for studying an interesting integrable per- 
turbation of the SW potential that has recently appeared in the literature known 
as the Tremblay-Turbiner- Winternitz (TTW) potential. Our interest will be in de- 
termining which potentials admitted by the TTW system remain multi-separable in 
the Euclidean plane. We remark that the results of the next two sections have been 
submitted for publication (see [57]). 
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3.2 Characterisation of the Smorodinsky-Winternitz 
Potential 

The superintegrability of the Smorodinsky-Winternitz (SW) potential is a simple con- 
sequence of its multi-separability with respect to both canonical polar and canonical 
Cartesian coordinates. The property of the coordinate web being canonical can be 
interpreted either geometrically, (see Figure 1.1), or with respect to the form of the 
general characteristic Killing tensor defining the orthogonal coordinate web (see [45]). 
Geometrically, canonical Cartesian coordinates are aligned with the coordinate axes 
and canonical polar coordinates have their singular point coincide with the origin of 
the coordinate axes; the associated Killing tensors for each canonical web respectively 
take the forms (3.6) and (3.8) given below. 

Since the SW potential is of the form V = f{x) + g(y), it immediately follows that 
the system enjoys separability with respect to (canonical) Cartesian coordinates, and 
so the Hamiltonian flow defined with respect to (3.4) admits a first integral quadratic 
in the momenta according to 



where K[ J is the Killing two-tensor whose normal eigenvectors generate the (canoni- 
cal) Cartesian coordinate web: 



readily integrate to obtain that Ui(x, %/) — -% — uj 2 x 2 . This provides us with the first 
non-trivial quadratic first integral. 

If we transform the potential of (3.4) to polar coordinates we immediately see that 
it is of the from V(r,6) = + g( r )i an d so that indeed the Hamiltonian system 
defined by (3.4) is also separable with respect to (canonical) polar coordinates. Thus, 



F 1 = iK{ i p i p j + U 1 , i,j = 1,2, 



(3.5) 




(3.6) 



Imposing the compatibility condition (3.3) requires that dlli = -ftTidV^, which we can 
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we have that it admits a second quadratic first integral according to 



F? = ±Ki j PiPj + U 2 , i,j = 1,2, 



(3.7) 



where K 2 is given by the Killing two-tensor whose normal eigenvectors generate the 
(canonical) polar web: 



we have obtained two non-trivial quadratic first integrals. 

The functional independence of F\ and F 2 follows from considering that the wedge 
product of their differentials with dH at a point p G E 2 is non- vanishing 



a criterion which yields that the set of Fi and F 2 defines a coordinate system, since 
this is equivalent to saying that the Jacobian determinant of the functions F\ and 
F 2 at a point p G M. is non-vanishing, and so by the Implicit Function theorem 
the set indeed forms a coordinate system (e.g. see [40]). Furthermore, since both 
are mutually in involution with the Hamiltonian, i.e. {F\,H} = {F 2 ,H} = 0, then 
we have at once that the Hamiltonian system (3.4) is superintegrable (since it is 
multi-separable) with 2n — 1 = 3 functionally independent first integrals of motion, 
namely the set {H, F\, F 2 }. Since the function F% = F\ + F 2 is also a first integral 
of motion, with an associated Killing two-tensor K\ + K 2 = Keh, where Keh 
denotes the (canonical) Killing tensor whose normal eigenvectors generate the elliptic- 
hyperbolic coordinate web, then it follows that the SW potential will also be separable 
in (canonical) elliptic-hyperbolic coordinates. Indeed, the canonical characteristic 
Killing tensor that generates the elliptic-hyperbolic coordinate web is given by 



(where k 2 is an invariant to be defined below) which makes it clear that any potential 




(3.8) 



In this case, we find upon integrating dU 2 = K 2 dV that U 2 (x, y) = + and so 



dH A dF 1 A dF 2 (p) ^0, 




(3.9) 
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of a Hamiltonian system compatible with both the valence two Killing tensor of polar 
type and the valence two Killing tensor of Cartesian type will also be compatible with 
the valence two Killing tensor of elliptic-hyperbolic type, by linearity of the exterior 
derivative. 

We can now make use of joint invariants, introduced in Section 2.4.3, to conclude 
that separation in (canonical) polar and (canonical) elliptic-hyperbolic coordinates 
implies that the quadrilateral (see Figure 2.1) S1S2S4S3 degenerates in the following 
manner. Denote by Si and S2 the singular points of the Killing two-tensor K3 whose 
first integral is given by F 3 as defined above (i.e. K 3 = Keh), an d similarly, take £3 
and £4 as the singular points of K 2 - Then since both Ki and K 2 , and hence K 3 
are in canonical form, we find that the quadrilateral S1S2S4S3 degenerates according 
to the invariant conditions S3 = £4, and dist(£i, S3) = dist(£2, S4). In terms of the 
fundamental joint invariants (2.61), we can equivalently write these conditions as 

A' 3 = 0, Ag = A 9 , 

respectively. 

Conversely, if we begin with a general potential V of a natural Hamiltonian 
(3.1) and impose that the Hamiltonian system be multi-separable with respect to 
canonical polar and canonical Cartesian coordinates, i.e. impose that it admits two 
quadratic first integrals given by Fi and F 2 as above, then integrating the condi- 
tions dUi = K{dV for i — 1,2 will yield the Smorodinsky-Winternitz potential, 
which follows from [22, 1]. Therefore, we can conclude that the superintegrable 
system given by the Smorodinsky-Winternitz potential is characterised by a pair 
of Killing tensors {K21K3) G /C 2 (E 2 ) x /C 2 (E 2 ) whose position in the orbit space 
(/C 2 (E 2 ) x /C 2 (E 2 ))/ 'SE(2) is determined by the invariant conditions 

A'^0, A' 3 = 0, A' 4 ^0, A' 6 ^0, A' 7 = A' 8 = A' 9 . 

Thus, we have the following: 
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Theorem 3.10. Let V be the potential of a natural Hamiltonian 

H(q\ Pi ) = \g ij W)PiPj + VtfUJ = 1,2 

satisfying d (^KgdV^j = 0, £ = 2, 3 for a pair of Killing tensors (K2, K3) G /C 2 (E 2 ) x 
/C 2 (E 2 ). Then the following statements are equivalent: 

(1) The pair (K2,K 3 ) G /C 2 (E 2 ) x /C 2 (E 2 ) is invariantly characterised by the con- 
ditions 

a; ^ 0, A3 = 0, a; ^ o, 

A' 6 ^0, A' 7 = A' 8 = A' 9 . (3.11) 

(2) The potential V given by the natural Hamiltonian is the Smorodinsky- Winternitz 
potential given in Cartesian coordinates (q l ,q 2 ) = (x,y) by 

V(x,y) = -u 2 (x 2 + y 2 ) + - 2 +P- 2 . (3.12) 

x z y z 

Thus, we have shown that by invariantly characterising the pair (K2, K3) G 
/C 2 (E 2 ) x /C 2 (E 2 ), we have in fact characterised the SW potential which they define. 
Next, we study what happens to the form of the SW potential in the case that we 
weaken the invariant conditions (3.11). 

3.2.1 Weakening the Conditions 

We shall first suppose that the Hamiltonian system (3.1) admits two quadratic first 
integrals of motion according to (3.5) and (3.7). Further, suppose that the Killing 
tensor Kp generates (non-canonical) polar coordinates, and that the Killing tensor 
Keh generates (non-canonical) elliptic-hyperbolic coordinates. Without loss of gen- 
erality, and to facilitate calculations, we will assume that K eh is in the canonical 
form, so that the elliptic-hyperbolic web is generated in the canonical position. Then 
the singular points of the two webs form a general triangle AS1S2S3 — see Figure 
3.1. Considering the fundamental joint invariants (2.61), this arrangement induces 
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Figure 3.1: Triangle whose vertices are given by the singular points of two Killing ten- 
sors, one non-degenerate and the other degenerate, (K eh, K p) 6 /C 2 (E 2 ) x /C 2 (E 2 ). 
In this case, Keh is in the canonical form. 

the following invariant conditions for the pair (Kp, K eh) £ /C 2 (E 2 ) x /C 2 (E 2 ): 

A' 6 ^0,A' 7 = A' 9 ^A' 8 . (3.13) 

We then have the following components for the respective Killing tensors as follows 
[45] 

(y-b) 2 (x-a)(y-b)\ 
i m u\ i ^2 ], a,beR, (3-14) 
[x — a)(y — o) [x — a) I 
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and 

K i £ H =( y2 + k2 ~ x 2 v ),k 2 eR, (3.15) 

V — xy x I 

where k 2 is the invariant given by (half) the distance between the foci of the elliptic- 
hyperbolic web (see Section 2.4.3). We then ask the question of how this setup affects 
the arbitrary parameters in the Smorodinsky-Winternitz potential (3.12). 

We must first determine the freedom we have in the arbitrary parameters afforded 
by the vector space of Killing tensors. Taking Fi (3.5) and F 2 (3.7) as the first 
integrals of the Hamiltonian system defined by (3.4), the Killing two-tensor obtained 
from taking a linear combination as 

K g = dK P + c 2 K EH + c 3 g, ci, c 2 , c 3 G K, 

must also satisfy the compatibility condition with the Smorodinsky-Winternitz po- 
tential V admitted by (3.4). However, this raises a concern since K g depends on 
six arbitrary parameters given by k 2 ,a,b,ci,C2,c 3 , and six also happens to be the 



dimension of the vector space of Killing tensors /C 2 (E 2 ). This implies that K g must 



in fact be given by the general formula 



K = ( ( 3 1 + 2p 4y + foy 2 )—Q — 

OX OXg 

+ fax- B 5 y- B e xy)—Q— (3.16) 
+ (/B 2 + 2fax + fax 2 )-^ © 

from which it follows by the compatibility condition (3.3) that the potential V is 
trivial: V = constant. Therefore, a non-trivial potential will be obtained if at least 
one of the six parameters vanishes. Considering the geometric implications afforded 
if we retain both Killing tensors Kp and Keh, the metric g, and k 2 , the only real 
freedom that remains is to switch off the parameters a and b. This presents us with 
the following cases: 

(1) G ^0, b^O, 

(2) a = 0, 6^0, 
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(3) a^O, 6 = 0, 

(4) a = 6 = 0. 

The first case imposed that the SW potential would assume a trivial form. We 
can introduce a new fundamental invariant (replacing A 9 in (2.61), since with respect 
to the two orthogonal coordinate webs A' g = A' 7 ) given by the area of ASiS 2 S 3 . In 
terms of the fundamental joint invariants (2.61) and the area of AS1S2S3, we then 
have the following invariant characterisation. 

(1) A; ^ 0, A' 3 = 0, A 4 ^ 0, A' 6 ^ 0, A' 7 ^ 0, A'g ^ 0, the area of ASiS 2 S 3 + °- 

(2) A; ^ 0, A' 3 = 0, A 4 ^ 0, A' 6 ^ 0, A' 7 - A' 8 = 0, the area of AS^Ss ^ 0. In 
this case the AS1S2S3 takes the form of an isosceles triangle. 

(3) A; ^ 0, A' 3 = 0,A' 4 ^ 0,A' 6 ^ 0,A' 7 ^ A' 8 , the area of AS 1 S 2 S 3 = 0. We 
remark that this is the case when the singular points all lie on the 

(4) A; ^ 0, A 3 = 0, A' 4 ^ 0, A'g ^ 0, A' 7 = A'g, area of AS X S 2 S Z = 

We can next characterise the arbitrary parameters a and b as joint invariants 
themselves by writing them in terms of the fundamental joint invariants (2.61). To 
this end, we consider the general Killing tensor whose components with respect to 
Cartesian coordinates (g 1 ,? 2 ) are given by (3.16) and suppose that f3§ 7^ 0, A = 0, 
where 

A = (/3 4 2 - Pi + - Pi)f + 4(/3 6 /3 3 + /3 4 /3 5 ) 2 . 
If we consider the Killing tensor Kp, we can compare it with (3.16) via 

K ij = £g ij + mKf j , (3.17) 

where Kf- = g^g^Kp 1 , with Kp 1 given by (3.14). Comparing components on both 
sides of (3.17) then amounts to identifying 

«! + 2a 4 y + a e y 2 = £ + m{y — b) 2 , 
«3 — — a 5 y — a^xy = —m(x — a)(y — b), 
a 2 + 2a$x + a^x 2 = £ + m(x — a) 2 , 
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which admits the following identification 

— «5 — «4 



(a, b) 



«6 «6 



Through a bit of algebraic manipulation, we can then obtain an expression for each 
parameter in terms of the joint invariants (2.61): 



(A' 4 ) 2 (A' 8 -A' 7 )-iA' 



2^6 



and 



\ 



A' 7 



2A' A 



Solving the PDEs that come from imposing the compatibility condition (3.3) for the 
Killing two-tensor defined by (3.14) and the SW potential (3.12), we arrive at the 
following equation 

9/T . a(b — y)a b(—a + x)B 
u 2 (bx - ay) + 1 + — = 0, 3.18 

which puts into evidence the following relationship 

(1) a 7^ 0, b 7^ yields the trivial potential V = 0, 

(2) a = 0, b 7^ yields the potential V — — , 

o; 

(3) a ^ 0,6 = yields the potential V = — , 

(4) a = 0, 6 = yields the general Smorodinsky-Winternitz potential (3.12). 

Henceforth, we arrive at the conclusion that the more geometric structure the pair 
(Kp, Keh) G K?{¥?) x /C 2 (E 2 ) has (in the sense of placing the polar coordinate web 
in restricted positions), the more general the Smorodinsky-Winternitz potential (3.12) 
becomes. In each of the cases described above the corresponding pair (Kp, K EH ) E 
/C 2 (E 2 ) x /C 2 (E 2 ) belongs to a different orbit of the group action SE(2) O /C 2 (E 2 ) x 
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/C 2 (E 2 ) (see Section 2.4.2, eq. 2.49). In fact, our above analysis shows that these 
corresponding orbits have been explicitly distinguished by means of joint invariants 
(2.61) of the action (2.49). Considering our work as an extension of the results 
obtained in [3], we see that we have demonstrated that the use of joint invariants 
to classify the orbits of (/C 2 (E 2 ) x K? {¥?))/ SE{2) leads to a classification of the 
corresponding superintegrable potentials defined by quadratic first integrals in E 2 . 
We remark on the further development of this perspective in Chapter 4. 

3.3 The Tremblay-Turbiner-Winternitz System 

We now consider the Tremblay-Turbiner-Winternitz (TTW) potential introduced in 
2009 [65]. It is defined in polar coordinates by the natural Hamiltonian 

H = V l + -p r + - 2 v 2 e - uj 2 t 2 + \ ( — ^- + - T ^- > ) , (3.19) 
r r 2 r z \cos z kv sin kO J 

which demonstrates that its potential, 

is of the from V(r, 9) = ^r-+g(r) and so is separable in (canonical) polar coordinates. 
By the results of Section 2.3, we have that the system admits a first integral of motion 
taking the form 

F(q\ Pi ) = \K^^)p iPj + U(q% i,j = 1,2, 

where K % i is the Killing two-tensor which generates (canonical) polar coordinates, 
i.e. (3.8) and U is determined via dU = KdV. We then seek to answer the following 
question posed in the introduction of this thesis: 

For which values of k is the TTW system multi-separable? 

In other words, we seek to determine for which values of k the TTW potential (3.20) 
admits an additional (non-trivial) functionally independent quadratic first integral of 
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motion. 

To allow for arbitrary k G K, we cannot transform the TTW system back into 
Cartesian coordinates. Indeed, such a transformation would require the use of trigono- 
metric identities in cos k9 and sin kO, which require k to in general be an integer. 
Thus we must carry out all our calculations in polar coordinates. To this end we 
shall first rewrite the formula defining a general Killing tensor K g G /C 2 (E 2 ) in polar 
coordinates. 

On this note, we make use of the result that on a manifold Ai of constant curva- 
ture, generalised Killing tensors are formed as sums of symmetrised tensor products 
of generalised Killing vectors on A4 (see [18]). This means that on E 2 , we have 

K g = A ij X i QX j +B i X i ®R + CR®R, i,j = 1,2, (3.21) 

for some tensor A % \ vector B l and constant C. If we take as our basis of Killing 
vectors on E 2 

n d sin 6* d 
X i = cos 9- 



and let 



dr r 89' 

. „ d cos 9 d .„ 

X 2 = sm9— + — , 3.22 

or r 69 

R- 9 



( Aij )= ( I' i 3 |. & = {Pt,fo), C = fr, (3.23) 

P3 P2 



then we find upon substituting (3.22) and (3.23) into (3.21) and taking the symmetric 
tensor products then 

K p = A lj Xi © Xj + B { Xi QR + CRQR 
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and we derive that the components of the general Killing tensor K p e /C 2 (E 2 ) defined 
with respect to polar coordinates are given by 

Kf = 2ft cos 9 sin 9 + ft cos 2 9 + ft sin 2 9, 

or p^or ■ a i 2 & cos 29 2ft cos g sing 2ft cos g sing 
iC ~ = 2p 4 cos 9 + 2ft sin 6 1 H 1 , (3.24) 

sr ry ry ry 

T^22 o 2ft cos 9 sin g ft sin 2 ft cos 2 9 2ft sing 2ft cos g 

-^P = A? ^ 1 o 1 n 1 ■ 

r ry Zi ry £1 ry£> ry ry 

Substituting Kp and the potential (3.20) into the compatibility condition, 

d (KpdV^j = 0, (3.25) 

we arrive at a complicated trigonometric equation in terms of the arbitrary parame- 
ters k, Ai, A 2 , ft, i = 1, ... 6 (see eq. (C.l) in Appendix C). To solve this equation, we 
expand it over the following set M of trigonometric functions 

M ={1, cos(g), sin(g), cos(2g), sin(2g), 

cos((l + k)9), sin((l + k)9), cos((l - k)9), sin((l - k)9), 

cos((l + 2k)9), sin((l + 2k)9), cos((l - 2k)9), sin((l - 2Jfe)0), 

cos((l + 3A;)g), sin((l + 3k)9), cos((l + 3Jfe)0), sin((l + 3A;)g), 

cos(2(l + k)9), sin(2(l + jfe)0), cos(2(l - k)9), sin(2(l - fe)g), 

cos(2(l + 2A;)g), sin(2(l + 2A;)g), cos(2(l - 2fc)g), sin(2(l - 2Jfc)0), 

cos(2(l + 3k)9), sin(2(l + 3A;)g), cos(2(l + 3A;)g), sin(2(l + 3A;)g)}, (3.26) 

which is the set obtained by expanding the trigonometric functions in (C.l), ignoring 
any trig identities. The set M is linearly dependent when the arguments of the 
trigonometric functions coincide (see Appendix D), a property which depends on the 
values assumed by the parameter k. By plotting the arguments as functions over k 
and determining the intersection points of these functions, we thus determine that 
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the set (3.26) is linearly dependent for the following values of k: 




14 



1 



1 3 2 

± — ,±-,±- 
16' 4' 3 



(3.27) 



313231321 
-,±-,± — ,±-,± — ,±-,± — ,±-,±-. 
8 3 10 7 14' 5 16 5 7 



In the case when k does not equal one of the above values, the set M is linearly 
independent and so we can set each of the coefficients to zero in the trigonometric 
equations that we obtain upon the substituting (3.21) and (3.20) into (3.25) and 
expanding over the set M. This argument leads to an overdetermined system of 
polynomial equations which we can then solve for the variables fc,/3j, and Xj,i = 
1, ... ,6, j = 1,2. Before completing this computation we shall reduce the number of 
polynomial equations with the following argument. 

Since the TTW potential (3.20) is separable in polar coordinates then if it were 
to admit another quadratic first integral of motion it would have been defined by a 
Killing two-tensor which is an element of the 5-dimensional vector subspace 



which is the vector space (3.21) with the "polar" term factored out. Substituting 
this Killing tensor K s into the compatibility condition (3.25) for the TTW potential 
(3.20), we arrive a system of equations for which it is possible to expand over the basis 
L = {1, cos 9, sin 9, cos 29, sin 29}. This results in the following necessary condition 
on the parameters of the Killing tensor K s : 



Since all parameters are assumed to be real, the above holds provided (3^ = /3s = 0. 
Therefore we conclude that if the potential (3.20) is separable in another orthogonal 
coordinate system, other than polar, it can only be the Cartesian system of coordi- 
nates. We can now repeat the same argument outlined above with a simpler Killing 
tensor. Namely, the Killing two-tensor that defines Cartesian coordinates in a general 



K s = A %3 Xi X 3 + B l X { R, i,j = 1,2, 



(3l + (3 2 5 = 0. 
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position [45], 




cos 2 (0 - <f>) 
§rsin(2(0 -</>)) 



§rsin(2(0-0)) 
r 2 sin 2 (6 l - 0) 



) 



(3.28) 



where denotes the angle by which the Cartesian coordinate system is rotated (recall 
the Cartesian web in Figure 1.2). Substituting (3.28) and (3.20) into the compatibility 
condition (3.25), we arrive at the following equation 

esc 2 k6 (6k cos 2(0 - 0) cot kO + (4 + 2k 2 - 3k 2 esc 2 k6) sin 2(6 - </>)) A 2 ^ 2g x 
+ sec 2 fc0 (-6A; cos 2(6* - 0) tan k6 
+ (4 + 2A; 2 - 3A; 2 sec 2 k$) sin 2(0 - 0)) Ai = 0. 

We can then expand this over a subset of iV C M given by 

N ={l,cos(20),sin(20), 

cos(2(l + k)9), sin(2(l + k)d), cos(2(l - sin(2(l - k)6), 

cos(2(l + 2ife)0), sin(2(l + 2k)9), cos(2(l - 2ife)0), sin(2(l - 2k)9), 

cos(2(l + 3k)6), sin(2(l + 3Jfc)0), cos(2(l + 3k)6), sin(2(l + 3A;)0)} (3.30) 

which we determine to be linearly dependent for the following values of k 



which are indeed a subset of the values obtained in (3.27). Expanding the equation 
(3.29) over the set N, and assuming it to be linearly independent (i.e. excluding the 
values of k given by (3.31)), we arrive at the following (reduced) system of polynomial 



k = ±1, ±2, ±2/3, ±1/2, ±2/5 



(3.31) 
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equations in the variables 0, k, Ai, A2: 



(-2 + k)(-l + k) (Xi-X 2 ) 


sin 


20 


= 0, 


-2 + fe(— 15 + 23k)) (Ai - A 2 ) 


sin 


20 


= 0, 


\~Z + A^lD + ZOK) ) [Ai — A 2y 


sin 


A 
Z(p 


- u, 


(1 + fcj(2 + fcj (Ai — A 2; 


sin 


OA 

Z(p 


= u, 


(,— 1 + — 1 + ZfcJ (Ai + A 2y 


sin 


OA 
Z(p 


- u, 


f 1 1 ZA 1 1 1 oM M 1 \ ^ 
[L + K)[L + IK) \A\ + A2j 


sin 


OA 
Z(p 


- u, 


( — 1 + ZK)(L + Ik) \A\ + A 2y 


sin 


A 

Z<p 


= u, 


(-2 + k)(-l + k) (Ai-Aa) 


cos 


20 


= 0, 


-2 + Jfe(-15 + 23Jfe)) (Ai - A 2 ) 


cos 


20 


= 0, 


(-2 + jfe(15 + 23ife))(Ai-Aa) 


cos 


20 


= 0, 


(1 + A;)(2 + fc) (Ai-A 2 ) 


cos 


20 


= 0, 


(-1 + k)(-l + 2k) (Ax + A 2 ) 


COS 


20 


= 0, 


(1 + fc)(l + 2A;) (A1 + A2) 


COS 


20 


= 0, 


(-1 + 2A;)(1 + 2A;) (A x + A 2 ) 


COS 


20 


= 0. 



There is no non-trivial solution for which the above system is satisfied. Considering 
the set of linearly dependent functions, we substitute into equation (3.29) with the 
special values of k given in (3.31) to yield the following possible cases 

(1) k = ±l 
Imposes = 0. 

(2) k = ±2 

Imposes {Ai = 0, = nn/2} or, {A 2 = 0, = nir/A}, n EN. 

(3) k = ±2/3, ±1/2, ±2/5, 
Imposes Ai = A 2 = 0. 
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Thus, we conclude that k = ±1 are the only values of k for which the full TTW 
potential (3.20) admits two quadratic first integrals of motion. For these values of k, 
the TTW system reduces to the Smorodinsky-Winternitz system discussed in Section 
3.2. Hence, we have proven the following 

Proposition 3.32. The general TTW potential given by (3.20) is a multi-separable 
superintegrable potential only when k — ±1. 

The case k = ±2 yields the well-known Calogero system on the line (see e.g. [10, 
14]). The conclusion of our above argument is that only when either Ai or A2 vanishes 
does this potential admit separation in general Cartesian coordinates according to a 
non-canonical Killing tensor Kq whose components are given by (3.28). In particular, 
the compatibility condition d(KcdV) = will be satisfied for the potential 

22, ^1 , ^2 



V(r,9) = -wV + — + 



cos 29 sin 29 

when either Ai = and <j) — nn/2, or A2 = 0, and (f) = rnr/4. 

Our results demonstrate that this approach to defining an infinite family of su- 
perintegrable systems appears to destroy the geometry of the associated orbit space 
(/C 2 (E 2 ) x K, 2 {¥?))/ SE{2). In other words, it seems unlikely for systems defined in 
this manner to admit two quadratic first integrals of motion for multiple values of 
the introduced parameter. In 2010 [50], this method was utilised by Post and Win- 
ternitz to construct another infinite family of superintegrable systems dependent on 
a parameter k, namely a potential that appeared as a second-order perturbation of 
the Kepler-Coulomb potential. They then established an equivalence between this 
new system and the TTW system through a Stackel transform, which maps Hamil- 
tonians to Hamiltonians, and hence developed a direct relation between the integrals 
of motion and trajectories of each system. Indeed, it would be of interest to pursue 
this equivalence through the invariant theory of Killing tensors which could possibly 
facilitate the search for these Stackel equivalent systems on a manifold. 
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Chapter 4 

Conclusions and Outlook 



We have invariantly characterised the Smorodinsky-Winternitz superintegrable po- 
tential modulo the action of the isometry group SE(2) in terms of joint invariants 
of the associated Killing tensors by taking them as points in the product space 
/C 2 (E 2 ) x /C 2 (E 2 ). We can extend this result to the problem of classification of arbi- 
trary superintegrable potentials defined in spaces of constant curvature by classifying 
first any associated Killing tensors that define first integrals of motion for the class 
of superintegrable systems in question (e.g., superintegrable systems that admit only 
quadratic first integrals of motion). These Killing tensors are identified as elements of 
the corresponding product spaces, whose factors are the corresponding vector spaces 
of Killing tensors representing the associated orbits in the induced orbit space under 
the action of the isometry group in such product spaces. 

Additionally, we have proven that the general TTW potential is a superintegrable 
system with two functionally independent quadratic first integrals of motion only for 
k — ±1, that is when it reduces to the Smorodinsky-Winternitz potential. If we 
allow one of the arbitrary parameters in the TTW potential to vanish, then a second 
multi-separable system is obtained when k — ±2. Therefore we have established that 
the general superintegrability of the TTW system for other k values must come from 
a (non-trivial) higher order (i.e. not quadratic) first integral. 

The results contained within this thesis motivate the complete classification of 
Hamiltonian systems with two degrees of freedom, which are superintegrable in the 
sense of being multi-separable. The introduction of an additional arbitrary parameter 
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into a superintegrable potential appears to be very destructive, in the sense that it 
destroys the property of multi-separability for all but one choice of the arbitrary 
parameter. This hints towards a characterisation on the Euclidean plane for all 
superintegrable systems obtained as perturbations of known superintegrable systems, 
and further, the classification of all multi-separable superintegrable systems defined 
on the Euclidean plane. 

Furthermore, we can extend our results to studying superintegrable systems de- 
fined on different two-dimensional spaces, such as on the Minkowski plane or the 
2D sphere. In the case of the former, the study becomes complicated by the fact 
that there are nine orthogonal coordinate webs and instead of a discrete set of singu- 
lar points, there are singular regions admitted by complex-valued eigenvalues of the 
associated characteristic Killing tensor. In the case of the 2D sphere, the study is 
simplified as a result of both the additional symmetry of the manifold and its com- 
pactness. Indeed, there are instead only two orthogonal coordinate webs defined on 
the sphere, and the set of singular points remains discrete. It is of immediate interest 
to the author to carry out the classification of all multi-separable systems defined on 
spaces of constant curvature with two degrees of freedom. 
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Appendix A 
Compact Formulae 



A.l The Action of the Isometry Group 

We provide compact formulas for the equations appearing in Section 2.4. Of course, 
these equations have been produced in the literature for higher (n > 3) dimensional 
spaces (see, e.g [29]), though we rewrite them specifically for E 2 at the indulgence of 
the author. 

This begins with the compact formulae for the group action SE(2) /C 2 (E 2 ), for 
the general Killing tensor given by 

K = A ij Xi Xj + B'Xi R + CRQ R, 

where we can define the parameters by 

a= (a M B= ( ft \ c = ftt 

\ A ft / V -ft / 

The transformation laws for the basis vectors are given by 

Xi = A j iXj, R = n j Xj + R, (A.l) 

where /i- 7 = e V'A^, where e k i is the two-dimensional Levi-Civita symbol with the 
usual orientation, e 1 2 = +l,e 2 i = —1, 5 l = (pi,p2) and the A J fc are elements of the 
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special orthogonal group, £0(2), i.e. 



A 



cos 9 — sin 9 



sin 9 cos 9 



) 



we find that the action of SE(2) /C 2 (E 2 ) induced by the action of SE(2) O E 2 , 
i.e. q % — > A 1 $ + 8*, can be written in a nice, compact way: 



where the parentheses around the indices denote symmetrisation as defined in Section 
1.3. This is equivalent to equation (2.49). 

A. 2 Kogan's Recursive Approach 

We can also provide compact formulae for those appearing in the derivation of group 
invariants using Kogan's recursive approach (see Section 2.4.2). Recall that we first 
compute the action of the subgroup of translations on the vector space of Killing 
tensors /C 2 (E 2 ). This is induced by the action of the translational group on E 2 , namely 
q l — > q l + 5 % , which yields the compact form of eq. (2.50) as 



The equivalent compact form of the normalisation equations (2.52) are then given by 



Finally, we needed to compute the action of the subgroup of rotations SO (2) O 
/C 2 (E 2 ). This action is obtained in a straightforward manner as the induced action 



A ij = AfcWpA^ + 2B*A fc V } + / W, 
C = /3 6 , 



(A.2) 



A lJ = A ij + 2B {i e ]) m 5 m + Ce i m e i ( .8 m 8\ 
B { = B j + Ce'tS', 
C = f3 6 . 



(A.3) 



<f = e % jB 3 /C. 
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of 50(2) E 2 , namely q l — > A^g- 7 to yield the following compact form of eq. 2.54 as 

A ij = A^A^A^, 

& = A/B l , (A.4) 
C = fa. 

This concludes the compact representation of the Section 2.4.2 equations. Though, 
we remark that the development can indeed be continued in order to derive the 
invariants of the vector space of valence two Killing tensors under the action of 
the isometry group, SE{2). Specifically, by substituting for 5 l into (A. 3) and using 
the transformation law given by (A.4). Provided that C ^ 0, then a basis for the 
invariants are admitted by the trace and determinant of A % K The case when C = 
must then be treated separately (see [29] for explicit details). 
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The Killing Tensor Equation in 
Local Coordinates 

With respect to local (Cartesian) coordinates, i.e. gij = 8{j, the Killing tensor equa- 
tion in Section 1.3, namely 

[g,K] = o, 

is equivalent to 1 

K (ij . k) = 0, i,j,k = 1,2, 

where the covariant derivative (see Section 1.3) is in fact a partial derivative since 
all the connection coefficients vanish in a space of constant curvature. Using the fact 
that Kij is symmetric, i.e. = Kji, and applying the symmetrization rule, this 
equation becomes 

Kij,k + Kj k>i + K ki j = 0, 
which yields a system of PDE's in K^, for i, j = 1,2. Namely, we have 

(1) K 11A = 0, 

(2) K 22 , 2 = 0, 

(3) K 12A = -1^11,2, 

lr This equivalence holds true only in E 2 , since contravariant and covariant components with 
respect to Cartesian coordinates are equivalent. 
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(4) tfi 2)2 = -§#22,1 

We will take the Killing tensor components as functions in the Cartesian coordinates 
(q 1 , q 2 ) = (x, y), so that, for instance #n,i = 9Kl Q^' y ^ ■ Integrating (1) with respect to 
x and (2) with respect to y then yields that Kn(x,y) = F(y) and #22(2;, y) = G(x), 
respectively, where F and G are arbitrary functions. Taking a partial derivative of 
(3) with respect to x, we find that 



8 2 K 



12 =0, 



dx 2 

ldF 

K 12tl (x,y) = H(y) = — . 



Integrating #12,1 with respect to y then yields that 

F(y) = -2 J H(y)dy + fa. (B.l) 
Analogously, if we take a partial derivative of (4) with respect to y , we find that 



d 2 K 



12 =0 



dy 2 

^# 12 , 2 = M(*) = ---, 



so 



G(x) = -2 J M(x)dx + (3 2 . (B.2) 

If we now use the fact that the partial derivatives commute, so f^J 2 = f^p, then 
since #12,1 (x,y) = H(y) and #12,2 = M(x), we find that 

OH _dM 
ay ox 

where f3 6 is a separation constant. Thus, we can separate these two PDEs to obtain 
that 

H{y) = -fay - (3 4 , M(x) = -fa - fa. 
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B The Killing Tensor Equation in Local Coordinates 

Substituting these into (B.l) and (B.2), respectively, we find after integrating that 
F(y) = fay 2 + 2 fay + fa, G(x) = fax 2 + 2f3 5 x + fa, 

so that 

K n = fa + 2/3 4 y + f3 6 y 2 , K 22 = f3 2 + 2(3 5 x + f3 6 x 2 . 

Lastly, integrating Ki 2 ,i(x, y) = H(y) with respect to x, and using -^12,2 = M(x), we 
will find that 

Ki 2 = /? 3 - PiX - (3 5 y + (3 6 xy. 

This gives the general solution to the Killing tensor equation defined on the Euclidean 
plane. Indeed, we see that it depends on 6 arbitrary parameters, fa, i = 1, . . . 6. 
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Appendix C 



A Complicated Trigonometric 
Equation 

Imposing the compatibility condition d (^KdV^j = in the case of a general Killing 
tensor K in six arbitrary parameters and the potential of the TTW system given in 
Section 3.3 as eq. (3.20), we arrive at the following trigonometric equation 

r V (r sin 20 (ft - ft) - 2r cos 20ft + 2 cos 0ft + 2 sin 0ft) 
+ 3k 2 r 4 sec 4 k9 (sin 29 (-ft + ft) + 2 cos 20ft + 2r cos 0ft 
+ 2r sin 0ft) A x + sec 2 fc0 ( (4 + r 2 + 2fc 2 r 4 ) sin 20 (ft - ft) 

- 2 (4 + r 2 + 2A; 2 r 4 ) cos 20ft 

- 2r (3 + 2A; 2 r 4 ) cos 0ft - 2r (3 + 2A; 2 r 4 ) sin 0ft) Ai 

+ 2A;r 2 (l + r 2 ) cot k9 esc 2 k9 (cos 20 (ft - ft) + 2 sin 20ft 

+ r sin 0ft - r cos 0ft) A 2 + 3A; 2 r 4 esc 4 k9 (sin 20 (-ft + ft) 
+ 2 cos 20ft + 2r cos 0ft + 2r sin 0ft) A 2 

+ -i + L 2 ^ (-< 4+r2+2fcV ) 5inM ' ft - fe » 

+ 2 (4 + r 2 + 2A; 2 r 4 ) cos 20ft + 2r (3 + 2A; 2 r 4 ) cos 0ft 
+ 2r (3 + 2£; 2 r 4 ) sin 0ft) A 2 -2A;r 2 (l+r 2 ) sec 2 k9 (cos 20 (ft 
+ 2 sin 20ft + r sin 0ft - r cos 0ft) Ai tan k9 = 0, 

which must be solved for co, X±, A 2 , k, ft, i = 1, . . . , 5. 



(C.l) 



-A) 
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Appendix D 



Linear Independence of a Set of 
Functions 



We provide a proof to the claim in Section 3.3 that the set M (3.26) is linearly 
independent when the arguments of the trigonometric functions do not coincide. 
Indeed, it suffices to prove the following. 

Proposition D.l (cf. Yuan [72]). The set of trigonometric functions defined by 



for positive r e M is linearly independent over K. 

Though it is straightforward to compute the Wronskian of the set we provide an 
alternative proof which allows us to extend the result. 

Proof. Suppose that we have a non-trivial dependency relation with a minimal num- 
ber of terms in sin(r6 l ) and cos(r#), 



Consider the largest possible real r for which d ro is non- vanishing. If we take a 
large, even number of derivatives eventually the d ro coefficient will dominate all other 
coefficients in the other cosine terms. Suppose we take a considerable amount of an 



{1, sin(r#), cos(r#)} 
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D Linear Independence of a Set of Functions 



even number of derivatives so that this is the case. If we then substitute with 6 = 0, 
we eliminate all sine terms and are left with a large number which cannot be equal to 
zero. Therefore, it must be the case that no cosine terms appear in the dependency 
relation. 

Similarly, we can consider the largest possible real r\ for which c n is non- vanishing. 
We will then take a substantial amount of an odd number of derivatives until the 
coefficient of the cos(r!#) term dominates all other coefficients in the other cosine 
terms (where these cosines come from differentiating sines). Substituting with 6 = 0, 
we again find ourselves with a large number which cannot be equal to zero. Thus, it 
follows that there can be no sine terms in the dependency relation. 

What remains is that 1 is the only function which can appear in the non-trivial 
dependency relation, which is clearly impossible. Thus, we conclude that no such 
dependency relation can exist and so the set must be linearly independent. □ 

Prompted by the discussion in [72], we arrive at a generalisation of the above 
proposition. Considering that 



it suffices to prove that the set of functions {e"' e } for r6l are linearly independent 
over BL Further, we can show that the set of functions {e z6 } for a z G C are linearly 
independent over C. 

Proposition D.2. The set {e zd } for z e C is linearly independent over C. 

Proof. Suppose we have a non-trivial dependency relation with a minimal number of 
terms given by 



as well. By our assumption, there must exist a non-zero z G C with a zo ^ 0, and so 





J> 2 e 2e = 0. 



Differentiating this equation with respect to 6 then yields that 



z ^e ze = 
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D Linear Independence of a Set of Functions 

subtracting this from the previous equation we find that 

^2 (z - z )a z e zt = 

holds as well. But this is a new non-trivial dependency relation in fewer terms, which 
contradicts our original assumption. Therefore, the set of functions {e z9 } is indeed 
linearly independent over C. □ 

We comment that the linear independence of the above functions follows from 
a more general property regarding eigenvectors of a linear operator, where those 
with distinct eigenvalues are linearly independent. Considering that the exponential, 
sine, and cosine functions are all essentially eigenfunctions of differentiation 1 , one can 
prove the result holds in this more general context. 



1 A non-zero function / is called an eigenfunction of a linear operator A iff it satisfies Af = A/, 
for some scalar A called the eigenvalue. 
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